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Abstract. Sphere packing problems have a rich history in both mathematics and physics; yet, 
relatively few analytical analyses of sphere packings exist, and answers to seemingly simple ques- 
tions arc unknown. Here, we present an analytical method for deriving all packings of n spheres in 
K 3 satisfying minimal rigidity constraints (> 3 contacts per sphere and > 3n — 6 total contacts). 
We derive such packings for n < 10, and provide a preliminary set of maximum contact packings 
for 10 < n < 20. The resultant set of packings has some striking features, among them: (i) all 
minimally-rigid packings for n < 9 have exactly 3n — 6 contacts, (ii) non-rigid packings satisfying 
minimal rigidity constraints arise for n > 9, (iii) the number of ground states (i.e. packings with 
the maximum number of contacts) oscillates with respect to n, (iv) for 10 < n < 20 there are only 
a small number of packings with the maximum number of contacts, and for 10 < n < 13 these are 
all commensurate with the HCP lattice. The general method presented here may have applica- 
tions to other related problems in mathematics, such as the Erdos repeated distance problem and 
Euclidean distance matrix completion problems. 



1. Introduction 

We consider all configurations of n identical impenetrable spheres in M 3 that maximize the 
number of contacts between the spheres. 

Our interest in this problem arose through the following question: can one direct the self- 
assembly of colloidal particles into any desired structure simply by imposing a spherically sym- 
metric binding specificity on those colloidal particles'^ Colloidal particles are small, (nanometer 
to micron sized) spherical particles in aqueous solution. In recent years a myriad of methods have 
been developed to control the binding of particles to each other, thus causing them to self-assemble 
into clusters [501 El M, M CGS GU • 

The interactions between the particles typically have a range much smaller than the particle 
size; the potential energy of such a cluster is then simply proportional to the number of contacts. 
Structures with maximum numbers of contacts are thus the minima of the potential energy, and 
correspond to what will form in thermodynamic equilibrium. 

Controlling which structures will form thus becomes a problem of controlling which structures 
correspond to energetic minima, or equivalently to contact maxima. Although minimal-energy 
clusters have been catalogued for many different potentials, e.g. the Lennard- Jones potential [21] . 
they had not previously been calculated for hard (impenetrable) spheres. Thus there was no 
detailed understanding of what structures could self-assemble in this colloidal system. 

The motivating question of how to control the self-assembly of a colloidal system could thus be 
broken up into 2 parts: (i) Determine what structures can self-assemble, (ii) Determine how to 
direct the self-assembly of the system such that only the desired structure(s) form. In this paper, 
we address (i). In [HJ[6], we address (ii). 

1 By 'spherically symmetric binding specificity' we are referring to a situation in which colloidal particles bind 
isotropically, but not all colloidal particles can bind to one another. For example, imagine assigning labels to the 
colloidal particles; A, B, C, and so on. Particles sharing common labels will be able to bind to one another, but 
colloidal particles that do not share a common label can not bind. 
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It should be noted that, in addition to self-assembly, the structures of small clusters of atoms 
or particles bear directly on problems central to materials science and condensed matter physics, 
including nucleation, glass formation, and the statistical mechanics of clusters [TH] H3, HI] . For 
instance, the first step towards understanding the thermodynamics of a particular cluster system 
is to calculate the ground states as a function of particle number n [T7j . 

The structures that can self-assemble in this colloidal system can be determined by solving the 
mathematical problem of which structures globally or locally maximize the number of contacts 
between n spheres in M 3 - i.e. all structures in which either (i) no additional contacts between 
spheres can exist, or (ii) a contact must first be broken in order to form an additional contact. We 
will refer to such structures as packings. We formulate our problem by focusing on enumerating 
only minimally-rigid sphere packings, which we define as packings with at least 3 contacts per 
particle and at least 3n — 6 total contacts. Minimal rigidity is necessary, but not sufficient, for a 
structure to be rigid. We previously detailed the ground states of these packings as well as some 
of their interesting features in [7]. Here, we present the results and method more completely. The 
method we introduce combines graph theory and geometry to analytically derive all minimally- 
rigid packings of n spheres. We perform this enumeration for n < 10 spheres. Due to the large 
number of packings that must be evaluated, this analytical method is implemented computation- 
ally, and near n = 10 we reach the method's computational limitations. Finding scalable methods 
for enumerating packings at higher n is a significant challenge for the future. 

Already by n = 10, a number of interesting features set in. For n < 9, all minimally-rigid 
packings have exactly 3n — 6 contacts. The first instance of a non-rigid sphere packing that 
satisfies minimal rigidity constraints occurs at n — 9, and more such non-rigid packings arise at 
n — 10. The first instance of packings with greater than 3n — 6 contacts occurs at n — 10. We 
discuss the geometrical manner in which these maximum contact packings arise, and conjecture 
that maximum contact packings for all n in this system must contain octahedra. We provide 
preliminary evidence for this maximum contact conjecture for n < 20, and we show that the 
putative maximum contact packings of 10 < n < 13 are commensurate with the hexagonal close 
packing (HCP) lattice, but that maximum contact packings of 14 < n < 20 are not. Furthermore, 
we show that the number of packings containing the maximum number of contacts is oscillatory 
with n, and we discuss the origins of these oscillations. 

The set of packings we enumerate includes, as a subset, structures previously observed and 
described in the literature: for example, it includes all minimal-second-moment clusters reported 
by [M], packings observed experimentally through capillary driven assembly of colloidal spheres 
[2U HI], as well as the Janus particle structures observed by Hong et al. [29J. 

This problem is closely related to several unsolved problems in mathematics, such as the Erdos 
unit distance problem (a.k.a. the Erdos repeated distance problem), Euclidean distance matrix 
and positive semidefinite matrix completion problems, and 3-dimensional graph rigidity. Thus the 
method and results presented here may have direct bearing on these problems. 

The organization of this paper is as follows: In the next section we outline our mathematical 
formulation of the problem and describe our methodology for finding all minimally-rigid packings 
of a fixed n. We combine graph theoretic enumeration of adjacency matrices with solving for their 
corresponding distance matrices. The elements of these distance matrices correspond to the rela- 
tive distances between spheres in 3-dimensional Euclidean space, and thus yield the packings that 
correspond to those adjacency matrices. Analytical methods for solving such adjacency matrices 
that scale efficiently with n do not exist. Sections [3] and [4] thus derive a method with improved 
scaling: Section tt\ derives geometrical rules that map patterns in adjacency matrices either (i) 
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to a configuration of spheres, in which case the adjacency matrix is solved for its corresponding 
distance matrix or matrices; or (ii) to an unrealizable configuration, in which case no real-valued 
embedding in 3-dimensional Euclidean space of that adjacency matrix exists in which spheres do 
not overlap. We show how these geometrical rules, combined with adjacency matrix enumeration, 
can lead to a complete set of minimally-rigid packings. Each time a new pattern is encountered, a 
new geometrical rule must be derived; thus this part of the method requires new derivations at each 
n, and does not scale efficiently^} In section |4| we derive a single geometrical rule (the triangular 
bipyramid rule) that can solve all iterative adjacency matrices. An iterative adjacency matrix is 
an n x n matrix in which all minimally-rigid m < n subgraphs also correspond to packings. This 
part of the method applies to any n, and thus scales efficiently for all n. Most adjacency matrices 
at small n are iterative; therefore, this greatly reduces the number of geometrical rules necessary 
to derive a complete set of packings. Section [5] describes the set of sphere packings we find from 
our study. We provide analytical formulas for packings up to n = 7, and the set of packings 
for n = 8, 9, 10 are included in the supplementary information [8]. Section [6] summarizes notable 
properties of the packings, including how the number of contacts changes with n, the emergence 
of minimally-rigid structures that are not rigid, and the emergence of maximum contact packings 
that are commensurate with lattice packings. Section [7] summarizes the main roadblocks towards 
obtaining results at higher n and contains several ideas and conjectures therein, and also discusses 
extensions to dimensions other than 3, as well as relevance to other problems of mathematical 
interest. Section [8] provides some concluding remarks. 

2. Mathematical Formulation 

We begin by presenting a two-step method for enumerating a complete set of sphere packings 
that satisfies minimal rigidity constraints. The set of all packings of n spheres is a subset of all 
possible configurations of n spheres. Thus, to enumerate a complete set of sphere packings we 
(i) use graph theory to enumerate all n sphere configurations, and (ii) determine which of those 
configurations correspond to sphere packings. The sphere packings we consider here correspond to 
maxima (local or global) in the number of contact^} Provided step (ii) is exact, this method will 
produce a complete set of packings. However, current analytical methods do not scale efficiently 
with n, and are therefore ill-suited for step (ii). Thus, in sections [3] and |i} we use basic geometry 
to derive an analytical method that can exactly determine which configurations correspond to 
sphere packings. Because the number of possible configurations grows exponentially with n, this 
analytical process must be implemented computationally. We focus our search to only those sphere 
packings satisfying minimal rigidity constraints (> 3 contacts per sphere, > 3n — 6 total contacts), 
because doing this guarantees that if a graph has an embedding in 3-dimensional Euclidean space, 
it corresponds to a sphere packing; whereas graphs that are not minimally-rigid can have 3- 
dimensional embeddings without corresponding to packings due to a continuous degree of freedom 
that allows for the formation of another contact. 

We also note that spheres can be thought of as points (Fig. [l^,b), where points correspond to 
the centers of spheres, and we measure the distance between spheres as the distance between their 
centers. Throughout this paper we will use the words sphere, point, and particle interchangeably. 



This method allowed us to enumerate packings up to n = 8, while standard methods in algebraic geometry 
(using the package SINGULAR) allowed us to enumerate up to n = 7. 

3 i.e. either the configuration of spheres can not form an extra contact (global maxima), or if an extra contact 
can be formed, it first requires the breaking of an existing contact (local maxima). 
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2.1. Graph Theory Produces the Set of Possible Packings. A configuration of n spheres 
can be described by an n x n adjacency matrix. A, detailing which spheres are in contact: A^ = 1 
if the i th and j th particles touch, and Aij = if they do not. A system of n spheres has n(n — l)/2 
interparticle distances; the 2 possibilities (touching or not touching) per distance thus leaves 
2n(n-i)/2 Qiiff eren t ways of arranging contacts amongst the distances. There are thus 2 n(jl ~ 1 ^ 2 
possible adjacency matrices, each of which potentially corresponds to a packing. 

Figure [I] shows a packing of 6 particles, both as a sphere packing (Fig. [T^i) and as points 
connected by line segments (Fig. [Tjo). The adjacency matrix corresponding to this packing is 
shown in figure ^Dp. The set of possible packings can be enumerated by considering all adjacency 
matrices. For n = 6, there are 2 15 = 32,768 different adjacency matrices. Table [T] shows that the 
number of adjacency matrices grows rapidly with n, reaching 3.5184 x 10 13 by n = 10; however, 
many of these correspond to the same structure due to particle labeling degeneracy. For example, 
switching labels 5 and 2 of figure [T] yields another adjacency matrix corresponding to the same 
structure^] Adjacency matrices corresponding to the same structure are isomorphic to one another 
- meaning there will exist a permutation of rows and columns that can translate one matrix into 
the otheiH 

To generate the complete set of possible packings, we need only enumerate nonisomorphic A's. 
Such algorithms exist; examples include nauty and the SAGE package called nice [121 EJ- The 
number of nonisomorphic matrices is much smaller but still grows exponentially with n. Table [T] 



shows this growth also; for example at n = 6 the number of potential structures is 156, and at 
n = 10 it is 12,005,168. 



Figure 1. Adjacency and Distance Matrix Representation of Packings. 

(a) 6 particle polytetrahedral sphere packing, (b) The same 6-particle packing shown 
in point/line representation, (c) Corresponding 6 particle adjacency matrix, A. (d) 
Corresponding relative distance matrix, T>. 

The set of *4's (potential packings) can be further reduced by imposing rigidity constraints. 
Most structures with less than 3n — 6 total contacts or less than 3 contacts per particle will not 
correspond to a packing because there will exist a continuous degree of freedom through which 
the structure can form one or more bonds. Rigidity requires (i) there be at least 3 contacts per 



Note that this is purely an exercise in switching particle labels, it is not a statement about the symmetry of the 
structure - we are not saying that particles 2 and 5 have the same contact distribution. The point is that many 
matrices can lead to the same structure, because how you label the particles is an arbitrary factor. 
5 See |27| for a nice example of such a permutation. 
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particle, and (ii) there be at least as many contacts as internal degrees of freedom - thus there 
must be at least 3n — 6 contacts^] 

Table [T] shows how imposing minimal rigidity constraints restricts the number of adjacency 
matrices. For n < 5 particles, this eliminates all but 1 adjacency matrix, thus identifying a 
unique packing for each of these n; the doublet, triangle, tetrahedron, and triangular bipyramid, 
respectively (section [5J. For n < 4, all relative distances within these packings are touching 
and are thus known. However, for n = 5, the packing contains one unknown relative distance, 
which must be determined. For n > 6, more than one minimally-rigid A exists, and thus rigidity 
constraints alone are insufficient to identify the number of sphere packings (or to solve for them, 
as more than 1 distance in each minimally-rigid A is unknown). 

2.2. Solving Potential Packings: Algebraic Formulation. To make further progress, we 
reformulate our problem algebraically. Each adjacency matrix element, Aij, is associated with an 
interparticle distance, 



yj ( x i ~ x j) 2 + (Vi - Vj) 2 + 0* - z j) 



2 



which is the distance between particles i and j whose centers are located at {xi } yi } z{) } {xj } yj } Zj) } 
respectively. The distances are constrained by the adjacency matrix as follows: 

(1) Aij = 1 Uj = 2r 

(2) Aij = rij > 2r 

where r is the sphere radius^} 

For adjacency matrices with 3n — 6 contacts, this leads to precisely as many equations as 
unknown^} The task is thus to solve for the > 2r given a particular set of = 2r. The 
particle configuration encoded by each A is thus specified by a distance matrix, V, whose elements 
T> — r 

The fundamental question is to find an efficient, exact method for mapping A — » T>. If any 
T>ij < 2r, this implies that particles i,j overlap; thus any T> with T>ij < 2r is unphysical. Figure 
[TJi shows the distance matrix corresponding to an n = 6 packing. The interparticle distance 
between each of the particles that are touching is normalized to 1; for this packing, this leaves 
three distances that need to be determined; ri 2 , r 23 , and r 34 . 

In solving A for V, the following scenarios are possible: 

(1) Continuous set(s) of real-valued V correspond to a given A, in which case the structure(s) 
are not rigid. 

(2) No real-valued V exists that solves A, in which case the structure is unphysical. 

(3) Finite, real-valued V exist that solve A. In this case, the structure(s) correspond to rigid 
sphere packing(s) provided that all T>ij > 2r. 



6 Note that in restricting to structures with exactly 3n — 6 contacts, we will also find structures with more than 
3rt — 6 contacts. This is because when we solve for the packings as outlined below, the solutions can end up having 
more contacts than are assumed in the algebraic formulation. 

7 Strictly speaking, Aij — implies that particles are not touching, and thus that rij > 2r. However, it is 
convenient to consider all solutions of A, which thus include the possibility = 2r for Aij = 0; in other words, a 
packing can be represented by multiple A's with different numbers of l's if the solution to Aij = forces it to a 1. 

8 Without loss of generality, we can set one particle to reside at the origin, another to reside along a single axis 
(for example the y-axis), and a third to reside in one plane (such as the xy-plane). 6 coordinates are then fixed, 
leaving 3n — 6 instead of 3n coordinates. 
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Table 1. The Growth of Adjacency Matrices with n. 



The number of adjacency matrices (constructed by [32]) decreases rapidly as isomor- 
phism and rigidity constraints are imposed. Iterative and non-iterative are defined 
in the text. The classification of whether an A is iterative or not is here shown 
after all rules for n — 1 particles are applied; thus the non-iterative column shows n- 
particle non-iterative structures only, and does not include non-iterative structures 
of less than n particles. Also note that the classification of whether an A is iterative 
or not is thus sensitive to which geometrical rules are used. The number of iterative 
and non-iterative As at n = 10 is different from [7] because some modifications were 
made to the code since the publication of that paper. Please see the supplemental 
information [8] for exactly which modifications have been made. 



For every nonisomorphic adjacency matrix, whether there exists a corresponding packing re- 
quires solving for T> and asking whether the resulting r^-'s satisfy these constraints. 

2.3. Limitations of Existing Solution Methods. The issue now becomes one of solving a 
system of n(n — 1)/2 constraints (3n — 6 equations and n(n— 1)/2 — (3n — 6) inequality constraints). 
Numerical approaches for solving such a system cannot be guaranteed to converge; for example, 
Newton's method requires an accurate initial guess for guaranteed convergence. When a solution 
does not converge, we do not know if it is because a solution does not exist, or because the initial 
guess is not sufficiently accurate^} Algebraic geometric methods (e.g. Grobner bases ) |TT] are 
effective, but these algorithms do not scale efficiently with n. Our own implementation^! of this 
method using the package SINGULAR [22] was only able to solve for structures up to n = 7. 

In the following section, we use another approach and derive a different geometrical method 
to efficiently solve for all sphere packings given a set of nonisomorphic, minimally-rigid A's. We 
implement this method up to n = 10, at which point we begin to hit some practical roadblocks; 
these are discussed at the end of the paper, where we outline potential ways to overcome them. 

2.3.1. Chiral Structures. Before proceeding further, it is worth remarking that structures with 
different handedness will correspond to the same distance matrix. We can analyze each V, and 

9 Furthermore, convergence to an unphysical solution lying within numerical error is also a problem. 

n = 8, using the software package SINGULAR [23], one matrix takes several hours to solve, and there are 
a total of 438 minimally-rigid non-isomorphic „4's. 
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determine whether it corresponds to a structure that has a chiral counterpart (see section 3.5). 
We refer to chiral structures as the same packing but different states - thus, a distance matrix 
having a left- and right-handed counterpart corresponds to 1 packing with 2 different states. Note 
that according to our definition, a different packing necessarily corresponds to a different state. 
Thus, the total number of states is equal to the total number of packings plus the total number 
of chiral counterparts. 

3. Geometrical rules solve for sphere packings 

We now show how geometrical rules can be used to effectively and analytically solve the class 
of polynomial equations that are generated by adjacency matrices. We use basic geometry to 
construct rules associating patterns of l's and O's in .A's with either a given relative distance, 
Dij, or an unphysical conformation (in which case no D > 2r exists). There thus exist two types 
of rules: Elimination rules eliminate an A as unphysical, and distance rules solve an Aij for its 
corresponding T>ij. 

3.1. Neighbor Spheres and Intersection Circles. With each sphere of radius r, we can as- 
sociate a neighbor sphere of radius R = 2r, whose surface defines where another sphere must lie 
if it is to touch the original sphere in question (Fig. [2^,). When 2 spheres touch, their neighbor 
spheres intersect in an intersection circle (Fig. [2)3). The radius of the intersection circle follows 
from straightforward geometry, and is ^R (see supplemental information for derivation [8]). 




Figure 2. The Neighbor Sphere. 

(a) A particle and its neighbor sphere. 4 particles are shown touching the center 
particle, and it is seen that their centers lie on the surface of the particle's neighbor 
sphere, (b) 2 touching particles. The associated neighbor spheres of the particles 
intersect in an intersection circle of radius (y/3/2)R. 



Now we can interpret each A^ = 1, in geometrical terms, as an intersection circle between 
spheres i and j. Minimal rigidity constraints imply that each particle is associated with at least 
3 intersection circles. Intersection circles can be used to derive geometrical rules because, in 
general, a packing of n particles involves intersections of intersection circles, and intersections of 
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intersection circles define points in space. A particle touching m other particles will lie at the 
intersection of those m neighbor spheres. For example, a particle touching the dimer depicted 
in figure will lie on the circumference of the associated intersection circle. The intersection of 
m > 3 neighbor spheres are points - and by defining points in space, basic trigonometry can be 
used to calculate the distances between those points, thus solving A's for Vs. 

3.2. Individual Geometrical Rules. Using intersection circles, we now derive several repre- 
sentative geometrical rules for eliminating and solving adjacency matrices. The supplementary 
information [S] contains the complete set of rules used to derive the results of this paper. 

3.2.1. Rule 1. The simplest rule arises from the fact that intersection circles can intersect in either 
2, 1 or points, but never in more than 2 points. This geometrical fact implies that any A with 
the following property is unphysical: any 2 of the set {Ajk, Aj P , Ak P } equal 1, and there exist more 
than 2 z's for which A%j = Aik = Ai P = 1. 

Physically, this implies that no more than 2 spheres can simultaneously touch 3 connected 
spheres. Three spheres are connected if at least 2 contacts exist between them (Fig. [3|. This 
in turn tells us how many identical spheres can mutually touch a trimer: 2. Figure [4] shows an 
example of an adjacency matrix that is unphysical for this reason: the blue highlighted section 
shows that spheres 4,5,6 make up a trimer; but the purple highlighted section shows that spheres 
1,2,3 all touch the same trimer. This is impossible given the argument outlined above and hence 
this adjacency matrix does not correspond to a packing. 




Figure 3. No More than 2 Particles can Touch 3 Connected Particles. 

Schematic of 3 linearly connected particles and their associated intersection circles. 
The center of each intersection circle lies at the midpoint of the line segment con- 
necting the associated points. There can never be more than 2 intersection points 
of these intersection circles, indicating that no more than 2 particles can touch the 
same 3 linearly connected particles. 
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Figure 4. Example of an Unphysical Adjacency Matrix, (a) An adjacency 
matrix that is unphysical because it implies more than 2 intersections of intersection 
circles. The blue highlights show that particles 4,5,6 make up a trimer. The purple 
highlighted part shows that particles 1, 2, and 3 all touch the same trimer, 4,5,6. 
(b) A sphere packing corresponding to this unphysical adjacency matrix (shown in 
both sphere and point/line representations). For it to be realized, 2 particles must 
occupy the same point in space. 



3.2.2. Rule 2. A trimer, a configuration of 3 spheres forming an equilateral triangle, is associated 
with 3 mutually intersecting intersection circles (Fig. [5^,). These 3 intersection circles intersect at 
2 points (shown in red). Here we calculate the distance between these 2 intersection points. 

Note that a particle lying at one of the intersection points forms the 4-particle packing (the 
tetrahedron). And that 2 particles, lying at each intersection point, form the 5-particle packing 
(the 5-point polytetrahedron). The distance between these 2 intersection points, h, is the only 
distance that is greater than R in the 5-particle packing (Fig. |6|. 

To calculate this distance, we note that the trimer and its associated intersection circles form 
the set of triangles shown in figure ^jp (where the dashed line indicates an out-of-plane triangle). 
We calculate a by considering the right triangle with sides \/3/2R — a, a, 1/2/2. Trigonometry 
then implies that a = R/(2y/3), and h = 2y2/3i£. 

This implies that the solution to an adjacency matrix corresponding to the 5-particle packing 

is 

/0 1 1 1 
10 11 
110 1 
1110 
\1 1 1 

where the right matrix is the corresponding distance matrix, V, and without loss of generality we 
have let R = 1. For n — 5, there is only 1 non-isomorphic minimally-rigid A. 
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We can formalize this construction as a distance rule, which can be used whenever a submatrix 
of some A has the same structure as the 5-particle packing. Such submatrices can be identified 
with the following pattern: Aij = Aik = Akj = 1, and there exist 2 points p for which A p i = 
A P j = A p k = 1. Whenever this pattern exists, the distance submatrix between the associated 
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(a) 




(b) 



Figure 5. The Intersections of 3 Intersection Circles. 

(a) A trimer and its corresponding intersection circles. The 3 intersection circles 
mutually intersect at 2 points, shown in red. (b) The triangles that relate the trimer 
to one of the points of intersection (red dot). The distance between this point and 
the center of the triangle is equivalent to half the distance between the 2 points of 
intersection (here denoted as h). 

points corresponds to T> for the 5-particle packing. In particular, this rule solves for the distance 
between the 2 points p, for example if p = I, m, then T>i m = 2-^/|i?. 

3.2.3. Rule 3. Another elimination rule follows directly from any distance rule, including Rule 2 
derived above. Suppose we determine that for a given pattern of Aij, the contact distribution 
implies that T>f. p > R. If it then happens that Ak P = 1, then this implies that all of the geometrical 
constraints cannot be satisfied simultaneously, so that A is unphysical. 

For example, if an A contained the intersection circle construction discussed in the previous 
section, that would imply that T>i m = 2a/2/3-R, but if the adjacency matrix also stated that 
Aim = 1, then that A would be unphysical. 

3.2.4. Rule 4- We can derive another set of geometrical rules by finding the maximum number of 
points that can lie on an intersection circle - this corresponds to the maximum number of spheres 
that can touch a dimer. Figure [7] shows the dimer (top and bottom spheres), as well as points 
lying on their intersection circle. 

The maximum number of spheres that can lie on an intersection circle is 5, and this can be 
calculated as follows: we divide the circumference of the entire intersection circle by the arc length 
swept out by 2 spheres lying a unit distance apartpj (see figure [7]). This arc length is given by 
S = rd, where r is the radius of the intersection circle, and 9 is the angle between 2 radial line 



Without loss of generality, we refer to the distance between two touching spheres, R, as the unit distance. 
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Figure 6. The 2 Intersection Points of a Trimer Correspond to the 5- 
Particle Packing. A 5-particle packing is shown with its point representation 
overlain. The center triangle of the point representation corresponds to a trimer, 
and the 2 points that contact the trimer correspond to the 2 intersection points of 
the trimer's 3 intersection circles. The 2 intersection points are shown in red, and 
h corresponds to the distance between them. 



segments. The law of cosines then implies that 
(3) 6 = cos" 1 ' 1 



3 / 

so that the number of points a distance R apart that can fit on an intersection circle is given by 

4 -= 2 - « 5.1043. 

fi?cos-(§) 

This indicates that (i) any A implying that more than 5 points lie on an intersection circle, and 
(ii) any A implying a unit distance between all m < 5 points lying on an intersection circle is 
unphysical. We can identify 5 points lying on an intersection circle by the following adjacency 
matrix pattern: Aij = 1, such that there are 5 points k for which Aik = Ajk = 1. 

To solve for the structure of m < 5 points lying on an intersection circle, we must compute 
the distances between the non-touching particles on the intersection circle (Fig. [71). Of these 
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Figure 7. 5 Points on an Intersection Circle. 

The intersection circle shown in black corresponds to the dimer in the center. 5 
points are shown lying on the intersection circle; this corresponds to 5 particles 
touching the center dimer. The radius of the intersection circle is (\/3/2)i2 (shown 
as dashed black lines), and connects the center of the dimer (which is the origin of 
the intersection circle) to points 1-5 on the intersection circle. The arc length swept 
out by one pair of particles on the intersection circle is S. It can be seen that the 
1st and 5th particles nearly touch. The space between them is not big enough to fit 
another particle, and thus it can be seen that no more than 5 particles can touch a 
dimer. 
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distances, we have already calculated that between points 1 and 3 and shown that it is = 2^2/3R 



(section 3.2.2). All of these distances can be obtained by the isosceles triangle with equivalent 
lengths y3/2R (corresponding to the dashed black lines in figure [7] - note that only 2 such lines 
are shown, but that they exist between the midpoint of the dimer and every point along the 
intersection circle). The unique length of the isosceles triangle will be the unknown distance, r^-, 
and the angle between the two y/3/2R sides connecting particles i and j will be called <pij. Thus, 
the unknown distances will all be given by 

(5) sin 
where 



1 \ -r- 
1 > 2 % o 



2 r V 

/ 2 



26 

2tt - 39 
2tt -49 



<Pl3 
014 
015 

where 9 is given by equation [3j thus yielding 

(6) r 13 = 2^-R 

(7) r 14 = 

(8) r 15 = ^R 

These calculations apply to any adjacency matrix in which Aij = 1, there exist n points, k, for 
which Aik = Ajk = 1, and there also exist n — 1 instances where A pq = 1 amongst the n points, 
k; where n = 3,4,5 for r 13 ,r 14 ,r 15 , respectively. Then the distance between the two endpoints of 
the n particles is given by V pk = r 13 , r 14 , r 15 , respectively. 
That is for 

p, q, I, then A pg = A g i = 1, and the distance V p i = 2^J^R. 
p, q, I, m, then A pq = A q i = A\ m = 1, and the distance V pm = |_R. 
p, q, I, m, z, then A pq = A q i = A im = A mz = 1, and V pz = ^R. 



n 


= 2: 


if k 


n 


= 3: 


if k 


n 


= 4: 


if k 



Note for n = 3, we have already identified this A pattern in section 3.2.2 whereas the n = 2, 4 
structures are new. Also note that, by symmetry, r 13 = r 24 = r 35 , and r 14 = r 2 5, and that these 
equivalences are identified by the above patterns in A. 

3.3. Using Geometrical Rules to Derive a Complete Set of Sphere Packings. The afore- 
mentioned geometrical rules can be used, along with the set of nonisomorphic adjacency matrices, 
to derive a complete set of packings for a given n. To explain more clearly how this is done, 
we show as an example the derivation for n = 7 particle packings. For this n, there exist 29 
minimally-rigid A's, all of which potentially correspond to packings (table [TJ. 

To these matrices, we apply the elimination rules just outlined, as well as those that appear in 
the supplementary information [S]. This immediately eliminates 24/29 ^4's as unphysical. Figure 
|8] shows which of the matrices are eliminated. Table [2] shows which rules are used to eliminate 
the .A's. 17/24 of the matrices are eliminated because they imply more than 2 intersections of 
intersection circles. Three of the matrices are eliminated because of the relative distance rule for 
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Figure 8. Eliminating Adjacency Matrices. There are 29 non-isomorphic 
adjacency matrices satisfying minimal rigidity constraints for 7 particles. 24 out of 
the 29 A's are eliminated by geometrical rules, which are shown here as color-coded 
X's - see table [2] and appendix A [8] for the corresponding rules. See figure [9] for 
solutions to the 5 physically realizable *4's (which correspond to the matrices that 
appear here without an X). 

three points on an intersection circle. Two matrices each are eliminated by rules acting on 5 rings 
and 4 rings, respectively. For the remaining five adjacency matrices, we apply distance rules to 
the A's to solve for the corresponding D's. Table [3] details which distance rules are used to solve 
for the packings corresponding to each A. The analytical solutions for the distance matrices as 
well as the associated packings are shown in figure |9j to each A (numbered by the order in which 
it appears in figure [8j in ascending order from left to right, and top to bottom, respectively), we 
apply the rules outlined in Table [3] to analytically solve for the packing. 

This is the set of 7 particle sphere packings. Note that the packing corresponding to graph 17 
(row 4 from the top) is the only one where distinct left and right handed structures are possible; 
thus it corresponds to 2 distinct states. 
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Figure 9. 7-Particle Packings. 

These are the 5 possible packings of 7 particles. Graph numbers appear to the left 
- corresponding to a sequential numbering of the graphs appearing in figure [8] from 
left to right and top to bottom, respectively. Following the graph numbers are A, 
V, and a picture of the corresponding packing, respectively. Note that graph 17 
is chiral - it has both a left-handed and a right-handed form (this can be seen by 
moving the top particle, currently shown on the left side, over to the right side). 
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Table 2. Elimination Rules Used for 7-Particle Packings. 

Each rule appears in its own section either in the text or in appendix A [8], where 
the complete set of rules is included. The rule column thus lists in what section(s) 
the relevant rules can be found. 



Color: 


Unphysical Because: 


Rule: 


X 


2 or more intersection circles intersect at more than 2 points 


A.l 

section 


3.2.1 


X 


All relative distances between 3 points lying on an 
intersection circle = R 


A.3 




sections 


3.2.2 3.2.3 and 3.2.4 


X 


A closed 5 ring surrounds a circle of intersection 


A.14 


X 


2 points on opposite sides of a closed 4 ring touch 


A.15 



Table 3. Rules Needed to Solve 7-Particle Packings. 

The rules listed here correspond to distance rules. Rule 'A.#' corresponds to rule 
# in appendix A [S], otherwise the relevant equation and section numbers are listed 
for rules found within the paper. (Note that rule 4, found in section 3.2.4 (eqn[6]), 
is the same as rules 2 (section [3.2.2 ) and A.l, rule 4 (eqn [7]) is the same as rule 
A. 2, and rule 4 (eqn[8]) is the same as rule A. 4). Graphs are numbered in ascending 
order from left to right and top to bottom as they appear in figure [8j These graphs 
correspond to the ones without X's. 



Graph Number: 


Rules Used: 


4 


section 


3.2.4 


eqns 


6. 


7 


and 8 


8 


section 


3.2.4 


eqns 


6 


and 7 


17 


section 


3.2.4 


eqns 


6. 


7 


and A. 11 


22 


A. 7 anc 


I A. 6 


26 


A.3 and A.9 



3.4. Non-Uniqueness of Geometrical Rules. Note that the geometrical rules described here 
are not unique in that (i) the rules themselves can be derived in different ways, and (ii) a different 
set of rules altogether could be derived/applied to solve for the same packings. This is simply one 
set of rules that works. One example of this is that either Rule 2 or Rule 4, equation [6] can be 
used to determine the unknown distance in a 5-particle packing. Another example of this is that 
all iterative packings can be solved using the triangular bipyramid rule that will be introduced in 
section [4] instead of using the aforementioned rules. There are undoubtedly many such examples, 
and the list of rules just presented were not derived with the goal of conciseness. They are 
complete in the sense that they allow one to solve all adjacency matrices for the n presented here. 
Beyond this however, they can only solve adjacency matrices containing the structure the rule 
identifies. If an A contains an identifiable structure as well as a non-identifiable (not previously 
encountered) structure, then the rules will solve the identifiable part and only partially solve that 
A for its corresponding T>. If an A contains no identifiable structure, then new rules must be 
derived to solve any part of it. It is because of these latter two cases that we continued to derive 
new geometrical rules as we increased in n. The triangular bipyramid rule in section [4] is a rule 
that we derived in order to have one general rule that could recognize a certain class of adjacency 
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matrix structures - the iterative class - for all n. Because it is a general rule, applicable to all n, 
its introduction makes the set of rules much more concise. 

3.5. Chirality. Once all packings have been derived by solving all *4.'s for their corresponding 
P's, we must determine how many states each packing has. If a packing is chiral, it will have more 
than one state. This will show up by a packing having a non-superimposable mirror image; for 
example a packing having different 'handedness,' such as distinct left and right-handed structures. 

One can calculate whether a packing is chiral as follows: The automorphism group of a packing, 
{a}, gives the set of self-isomorphisms - i.e. all possible permutations of the structure into itself. 
Each element of the automorphism group will thus correspond either to a rotation or to a reflection. 
Rotations are transformations with determinant = 1, and reflections are transformations with 
determinant = —1. Thus, one can construct the set of all isomorphic graphs, {£>}, and construct 
the automorphism group for any one V>i within the isomorphic set. If 3 any matrix Dj G {T>} that 
is isomorphic to T>i only through reflections and not through rotations - i.e. if the isomorphism 
group of T>j to T>i corresponds to transformations that all have determinant -1, then the packing 
T> is chiral. 

A property related to chirality is the symmetry number, a. This corresponds to the number of 
ways that a structure can be rotated into itself. The symmetry number is necessary for calculating 
the equilibrium probability distribution of packings S3] . The symmetry number of a packing, 
T>, will be equal to the number of transformations within the automorphism group of T> that have 
determinant 1. Thus, if a packing has no reflections, then the symmetry number will simply equal 
the size of the automorphism group. If a packing has reflections, then the symmetry number will 
equal the size of the automorphism group divided by 2 (dividing by 2 will remove all automorphism 
mappings corresponding to reflections, and not rotations). 

Related to both symmetry numbers and chirality are point groups. A point group is a group 
of symmetry operations which all leave at least one point unmoved. Point groups have been 
calculated for many structures - and there exist programs that allow one to enter in a set of 
coordinates and retrieve the point group corresponding to those coordinates |46j. Symmetry 
numbers and chirality can alternatively be calculated directly from the point group of a structure. 
For example, compounds in the C m point group, where C m is the cyclic group consisting of 
rotations by 360°/m and all integer multiples (where m is an integer), are always chiral [25] . 

Point groups, symmetry numbers, and chirality of packings are included in the lists of packings 
appearing in section [5] and in the supplementary information [8] . The growth of chiral structures 
with n is interesting - surprisingly, over half of all 9-particle packings are chiral - see table |4j 

4. One Geometrical Rule That Solves for all Iterative Packings: The 

Triangular Bipyramid Rule 

In principle, these types of geometrical rules can be used to derive a complete set of sphere 
packings for any number of particles, n. However, in practice, the number of rules used here 
grows too quickly with n for this to be a practical method: at n = 5 spheres, only 1 rule is 
required, 3 rules are needed at n = 6, 12 rules at n = 7, and at n = 8, 14 rules solve 435/438 
minimally-rigid non-isomorphic A's. This leaves 3 unsolved *4.'s for which more geometrical rules 
must be derived; looking ahead at the 13,828 and 750,352 A's that must be solved at n = 9, 10, 
respectively, it becomes clear that deriving a rule or set of rules that does not grow significantly 
with n is a necessary step. Here, we derive one geometrical rule that can solve one class of packings 
for any n, thereby greatly reducing the number of rules needed to derive a complete set of n sphere 
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(a) (b) 

Figure 10. Iterative Packings. Two examples of iterative packings, (a) A 6 
particle polytetrahedron (red) with one particle added to it (blue). This decomposes 
into a tetrahedron (blue) plus a 6 particle polytetrahedron (red), with a shared 
triangular base (purple), (b) 2 joined octahedra (one red and one blue, with a 
shared purple triangular base) forming a 9-particle packing. 



packings. In section 7.2 we discuss how one geometrical rule can also be used to solve the other 
class of packings for all n. 

Packings can be broken up into two types or classes: iterative and non-iterative, or new seeds. 
Iterative packings are n-particle packings that are solely combinations of packings of less than n 
particles (see Fig. [lO|. New seeds are n-particle packings that cannot be constructed solely out of 
packings of less than n particles - i.e. they contain within them (in part or in whole) an inherently 
new structure (Fig. 11). Put another way, iterative packings correspond to *4's for which all 
minimally-rigid m x m subgraphs, m < n, correspond to packings that have been identified at 
lower n. 
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Figure 11. New Seeds. The octahedron is an example of a new seed, (a) The 
'base' of an octahedron has a continuous degree of freedom through which a 5 particle 
polytetrahedron can form, and thus is not a packing. The continuous degrees of 
freedom are shown as dashed lines; bringing either of the pairs of particles connected 
by these dashed lines into contact forms the 5 particle polytetrahedron. (Note that 
the 5 particle structure shown in (a) is not minimally-rigid as it has fewer than 
3n — 6 = 9 contacts.) (b) Once a 6th particle is added, the octahedral structure can 
be stabilized, thereby forming a new seed. This 'new seed' is an inherently new 
structure. 
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4.1. Solving iterative structures. An iterative packing is a polyhedron containing solely pack- 
ings of less than n particle^} Thus any iterative packing can be decomposed into 2 joined 
polyhedra (see Fig. [To] - the red and blue packings are the joined polyhedra). The 2 polyhedra 
are joined via a common base of particles (shown in purple )p} Because the joined polyhedra are 
less than n-particle packings, all of their mtrapolyhedral distances are known from lower-order 
packings. 

Thus, deriving one geometrical rule that can solve for all iterative packings requires solving 
the following geometrical problem: Given 2 joined polyhedra, where all znirapolyhedral distances 
are known, derive a general formula for the mterpolyhedral distances. Note that the solution to 
this problem immediately extends to unphysical iterative structures as well, as they are composed 
of structures of less than n particles, where either (i) one or more of the joined structures is 
unphysical, or (ii) the particular combination of the structures is unphysical. 

The geometrical problem is solved with the following observation: An explicit analytical formula 
for the distance between any 2 points can be derived if those 2 points can be related to a common 
triangular base. Let there exist two particles i,j whose interparticle distance, r^, is unknown. 
If there also exist 3 particles, k,p,q, with known interparticle distances (r kp , r kg , r pg ), and if the 
distances between i, j and the 3 particle base (r ip ,r ik ,r iq ,rj p ,rj k ,rj q ) are also known; then there 
exists an analytical relationship for the resulting r^. We call this the triangular bipyramid rule 
because the 5 points i,j, k,p, q together form a (potentially irregular) ditetrahedron or triangular 



bipyramid (see figure 19). We show the rule here, while a complete derivation can be found in the 



supplemental information (Appendix A [8]). 



The dihedral angles A2 and A 3 of figure [19] are given by: 

' cos a 3 — cos 63 cos c 3 



A 3 = cos 
Ao = cos" 



sin c 3 sin b 3 
cos 0,2 — cos 62 cos C2 



sin C2 sin 6 2 

These formulas are essentially obtained using the spherical rule of cosines - see Appendix A for 
details [8]. 

The dihedral angle A\ will be either the sum or the difference of the dihedral angles A 2 and 



A 3 (see Fig. 19), depending on whether the points i, j lie on the same or on opposite sides of the 
base p, k, q. If i, j lie on the same side, then Ai = \A 2 — A 3 \, and if i, j lie on opposite sides then 
A x = A 2 + A 3 . 

The angle a x is then given by 

(9) a\ = cos _1 (sinci sinfei cosAi + cosfei cosci) 

and from the law of cosines, we can then calculate r^ : 



(10) nj = r\ + r 2 pj - 2r ip r pj cos a x 

Associated with each we have 2 possible A x , and thus 2 possible solutions (similar, in 
principle, to one having 2 possible solutions to a quadratic equation). 

12 An iterative A is an n x n graph composed solely of to x m (to < n) subgraphs, each of which correspond to 
minimally-rigid ^4's of less than n spheres. 

13 Given the minimal rigidity constraints we have imposed, this common base will always consist of at least 3 
particles. 
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Figure 12. The Triangular Bipyramid. 

The triangular bipyramid (or ditetrahedron) constructed in the triangular bipyramid 
rule. The center triangle (kpq), shown in red, corresponds to the common 3 particle 
base. Particles % and j are related to one another through the common base. The 
distance between i and j, r^, shown as the dash-dot blue line, is unknown. a\ 
corresponds to ^jpi. A\ is the dihedral angle between Aipk and Ajpk, Ai is the 
dihedral angle between Aipk and Akpq, and A3 is the dihedral angle between Akpq 
and Ajpk. Points i and j can either both lie on the same side of the base kpq or 
each lie on opposite sides of the base (indicated by the dashed lines, that can either 
go into or come out of the plane). If i,j lie on the same side, then A\ is equal to 
the difference of A 2 and A 3 , and if i, j lie on opposite sides of the base then Ai is 
equal to the sum of A2 and A3. When all distances other than are known, then 
an explicit analytical formula can be derived to solve r^. 



4.2. Applying the triangular bipyramid rule. The triangular bipyramid rule can be used to 
solve all iterative packings as follows. We first search for subgraphs of A corresponding to lower-n 
seeds. The elements of T> corresponding to these lower-order structures are known and inserted 
appropriately. If A is iterative, all minimally-rigid subgraphs of m < n particles (i.e. m subgraphs 
with at least 3m — 6 contacts and at least 3 contacts per particle) will correspond to m-particle 
packings. Once all lower-order seeds are inserted as appropriate, all unknown correspond to 
the distances between the spheres of different known lower-order subpackings. The triangular 
bipyramid rule is then applied to each unknown elementp] of T>. For each unknown distance, ry, 
both solutions are potentially stored, as are all possible sets of unknown distances {7^}. Along 

14 If A is not iterative, there will exist unknown r,-j that do not correspond to distances between spheres of known 
subpackings. In this case, equation [32] will contain at least 1 unknown element on the right hand side and can not 
be applied directly. 
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the way, each solution is tested for consistency, and it is always possible that both, 1, or neither 
solution will be consistent. Once all locally consistent are stored, the resultant {%} are tested 
for global consistency. 

A solution will be inconsistent, and thus unphysical, for one of the following reasons: 

1. It violates the triangle inequality (meaning that no real solution exists - this shows up as 
the absolute value of the argument of the inverse cosine being greater than 1). 

2. One or more distance(s) are less than R. 

3. Different triangular bases lead to different r^; this indicates that a structure is in conflict 
with itself. One part of it implies it should have one structure, whereas another part implies 
a different structure. Such a structure is physically and mathematically inconsistent. 

Violation |2] arises within the 5 particles of one triangular bipyramid, and thus registers a physical 
local inconsistency. Violation [T] occurs within individual triangular bipyramids, as each is de- 
termined, in which case it registers a local inconsistency; as well as over the entire set of triangular 
bipyramids, once all {r^} have been determined, in which case it registers a global inconsistencjj^j 
Violation [3] occurs when solutions are consistent within individual triangular bipyramids, and thus 
locally consistent, but inconsistent within combinations of triangular bipyramids - these solutions 
are thus globally inconsistent. This violation can be checked as follows: Figure [19] shows that the 
dihedral angle, A\, is given by either the sum or the difference of A2 and A3, if particles i and 
j lie on opposite sides or on the same side of the triangular base, respectively (to within a 2ir 
modulation of course). Test all possible 5 particle combinations of triangular bipyramids within 
the n particle structure, and if there exists a triangular bipyramid that does not satisfy 

A 2 + A 3 

rm a - \ A i-M 

(U) Al ~ 2n-(A 2 + A 3 ) 

2tt- \A 2 -A 3 \ 



the solution is globally inconsistent. (In calculating in equation 32, we need not consider the 
latter two A\ solutions, as cos(2-7r — x) = cos(x)). 

There is one scenario in which 5 points need not satisfy this global consistency check, and that 
is when 3 or more of the 5 points lie in a line. In this case, the 3 or more points define a line and 
not a plane, and thus the dihedral angle is not defined, and equation [TT] can not be applied. This 
situation is encountered in certain structures that contain octahedra (for example, see graphs 5416 
and 10664 at n = 9 in the supplementary information |8J). In this scenario, the following global 
check can be performed: for m points lying in a line, there must exist the following number of 



15 In this case, some triangular bipyramids are locally consistent, whereas others are not. All possible triangular 
bipyramids of a structure need not be tested to solve for all r^-, thus it is important to check all bipyramids to 
ensure global consistency once {r^} have been determined. This violation is related to violation[3j except that here 
the violation is registered between the angles associated with the line segments, and in violation [3] the violation 
occurs within the dihedral angles. 
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line segments with a given minimum distance: 

Number of Line Segments Minimum Length of Line Segment 
m — 1 R 
m-2 2R 

(12) 

m — (m — 1) (m — 

Note that the first constraint, having m — 1 line segments with a minimum length of R, is auto- 
matically satisfied by the fact that we eliminate all solutions containing a distance < R. Up to 
n = 10, it turns out that performing this check for m = 3 alone is sufficient to find all globally 
inconsistent solutions. We evaluate that 3 points lie in a line by identifying that the angle associ- 
ated with those 3 points is or tt. We then check that there is at least one distance amongst the 
3 points that is > 2 (where once again, R = 1 without loss of generality). If this is not satisfied, 
then the packing is globally inconsistent and is eliminated. Note that the fact that some dihedral 
angles may not be defined was overlooked in [7], and is the primary source of the difference in 
numbers reported here in tables [I] and [4] versus the numbers previously reported in [7j. 

4.3. The Growth of New Seeds. For new seeds, we have a structure that contains an inherently 
new polyhedron, and thus some or all mtrapolyhedral distances are also unknown, i.e. one or 
more of the 9 distances within {ri P ,rik,r iq ,rj p ,rjk,rjq,r p k,r pq ,rk q } are unknown. Thus, deriving 
the equation for r^-, as is done for iterative packings, will yield one equation with more than one 
unknown. The triangular bipyramid rule, therefore, can not be applied directly to new seeds, and 
new geometrical rules must be derived to analytically solve non-iterative .4,'s. 

Using a general rule to solve for iterative *4's, and deriving individual geometrical rules for 
non-iterative A's is feasible so long as the non-iterative .A's do not grow too quickly with n. 
Table [I] shows that this is the case for n < 9, where the number of non-iterative .A's is 5 or less. 
However, at n = 10, there are 94 non-iterative ^4's^j To sift through the 94 potential seeds at 
n — 10 requires inventing new geometrical rules, and the growth of such rules demonstrates that 



the method we have described does not scale efficiently with n. In section 7.2 we will discuss a 
potential extension of the triangular bipyramid rule, which might be able to break this bottleneck, 
at least computationally. Here, we present the packing results derived from a combination of the 
triangular bipyramid rule and individual geometrical rules. For n < 9, we have analytically 
solved for all packings. At n = 10, we analytically solve for all iterative packings, and produce 
a preliminary list of new seeds, found by solving the non-iterative A's numerically using Newton 
Iterations^ 7 ! 



16 Note that the non-iterative and iterative „4's listed in this table are constructed after applying the geometrical 
rules for n < 8 that appear in the text and in the supplementary information [8]. Thus, this reflects the number 
of iterative and non-iterative .A's with respect to these geometrical rules, and not the absolute number of iterative 
and non-iterative ,/Ts 

17 Our runs of Newton Iterations used random initial guesses (between -5 and 5 for the coordinates of the 
particles). We performed 20 iterations of each initial guess, and approximately 150,000 total runs. Every matrix 
for which a solution was found, was found multiple times. We believe this to be a reasonably thorough search and 
that this preliminary list of new seeds at n = 10 might be complete. It is worth noting that the preliminary list of 
new seeds reported here found by Newton Iterations is the same as the list reported in [7] found by constructing 
the non-iterative A's manually with the construction toy Geomags. 
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5. The Set of Sphere Packings 

Here we present the list of sphere packings derived by this method^ 5 } In principle, the analytical 
method outlined here will yield a complete set of minimally-rigid packings. However, we have not 
implemented the triangular bipyramid rule symbolically, leading to the following practical issues, 
which could lead to numerical errors: 

Numerical round-off error: All calculations involving the triangular bipyramid rule are subject 
to numerical precision. Our algorithm eliminates many packings by finding situations where the 
argument of the cos -1 term is larger than unity; this can also occur erroneously due to round-off 
error, causing packings to erroneously be recognized as unphysical. Similarly, round-off errors 
are possible when checking for the consistency of a packing; in checking the equivalence between 
dihedral angles, or checking that there exists at least one distance greater than or equal to 2 when 
3 points lie in a line. Round-off issues could be improved by using general precision libraries such 
as grap and mpfr [H [2] , or altogether avoided by doing all calculations symbolically. Thus, while 
the analytical method presented here should in principle yield a complete set of sphere packings, 
practical issues such as these are a source of error. 

We present a complete set of sphere packings of n < 9, save round-off error. At n = 10, we 
present a complete set of iterative packings and a preliminary list of new seeds. Packings of n < 7 
particles are included here, and packings of n < 8 < 10 particles are included in supplementary 
information [8]. All results are summarized in table |4j 

In the list presented here, corresponds to the point group, and a to the symmetry number. 
We have included the 2nd moment of each packing, and a '*' appears next to the 2nd moment that 
corresponds to the minimum of the 2nd moment of n particles. The 'special properties' column 
denotes whether a structure is convex, a new seed, chiral, or non-rigid. If the special properties 
column is blank, then that packing contains none of these properties. 



18 The number of packings presented here for n = 9, 10 differs from the number we reported in [7]. This is 
primarily because our previous code did not run a check to ensure that the dihedral angle was well defined when 
checking for global consistency. As mentioned in section [472| this can occur when the 3 points used to define a plane 
are collinear, occurring in some of the packings that contain octahedra. In our original code (used for [7]), the 
dihedral angle check was still being performed in such an instance and erroneously deemed some packings globally 
inconsistent. It was a personal communication with Rob Hoy, who extended this work in (30] . that brought it to 
our attention that 2 packings were missed in our n = 9 list. In examining this discrepancy, we discovered that this 
issue with the dihedral angles was what caused these packings to be missed, and we have since made the relevant 
correction in the code. This caused 52 packings to be realized at n — 9 and many more packings to be realized 
at n = 10. We also made 2 more modifications to the code, such as further correcting for numerical round-off 
error, which further caused several more packings to be registered at n — 10. Please see supplemental information, 
Appendix C [8] for a complete list of the changes that were made to the code. 
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2 Particle Packings 



Packing 1 (Graph 1): 
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3 Particle Packings 



Packing 1 (Graph 1): 
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4 Particle Packings 



Packing 1 (Graph 1): 
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5 Particle Packings 

Packing 1 (Graph 1): 
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6 Particle Packings 

Packing 1 (Graph 2): 



.4 : 



/O 1 1 
OOOl 
10 
110 



V : 



v ■ 

( 

1 
1 
1 

V i 



1 1 



2Ji 



11110 1 



oj 



( 4/(3^3) \ 
1/3 

(2/3)v^73 


-1 


20/(973) 
-4/9 
(10/9) 




1/(2^3) 
-1/2 
^273 
V3/2 
-1/2 




) 




2nd Moment 





a 


Special Properties 


3.35185i?" 


Civ 


2 





Packing 2 (Graph 4): 
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7 Particle Packings 



Packing 1 (Graph 4): 
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Packing 2 (Graph 8): 
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Packing 3 (Graph 17): 
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Packing 4 (Graph 22): 
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Packing 5 (Graph 26): 
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6. Properties of Packings 

Here we highlight some interesting properties of packings. 

6.1. New Seeds. New seeds are interesting because they are inherently new structures of n par- 
ticles. They are also 'generating sets,' i.e. once they exist at a given starting n = m, they are 
propagated iteratively for all n > m. Geometrically, new seeds are inherently global structures, 
stabilized exactly by the n particles for which that new seed arises. Iterative packings are geomet- 
rically locally stable, in that subsets of less than n particles within the packing also correspond 



to packings. Thus, new seeds are unique events for a given number of particles, n. Figure [13 
shows all new seeds of n < 10 particles (where the set of new seeds at n = 10 is putative). The 
proportion of new seeds to total packings is relatively small for small n, which can be seen in table 

m 

6.2. Rigidity. We have enumerated packings satisfying minimal rigidity constraints; these con- 
straints are necessary but not sufficient for rigidity, and thus we can find non-rigid packings that 
satisfy these constraints. For these packings, there exists a degree of freedom in which particles can 
move without breaking or forming additional contacts. The first instance of a non-rigid packing 
occurs at n — 9, at which there is one. At n = 10, there are 4 non-rigid packings: 1 non-rigid new 



seed, and 3 iterative non-rigid packings that derive from the n = 9 non-rigid new seed (Fig. 14 



and 15) 



These 10 particle non-rigid packings will iteratively produce at least m > 1 non-rigid packings 
at n = 11, and so on. All non-rigid packings enumerated thus far contain at least 2 deformable 
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Table 4. Packings. 



Total number of packings found by the current study, compared to those found by 
Hoare and Mclnnes [27], who used only the tetrahedral (n = 4) and octahedral 
(n = 8) seeds to iteratively calculate hard sphere packings by adding one particle 
at a time to (n — l)-particle packings. They include chiral structures within their 
list of packings, so that a left and right-handed structure are considered to be 2 
packings. We distinguish between chiral structures and packings, such that a left 
and right-handed structure is considered to be 1 packing with 2 distinct states (for 
n < 10, left /right-handedness is the only type of chiral packing encountered). The 
number of packings having chiral counterparts is included in the column marked 
'chiral.' The total number of states per n is equal to the number of packings plus 
the number of chiral structures. This is included in the table, along with the number 
of packings corresponding to new seeds and to non-rigid structures. The number 
of packings and total states for n — 9, 10 differ from the numbers we previously 
reported in [7J primarily because in [7J we did not run a check to confirm that the 
dihedral angle was defined. Here we report 2 more packings (and 3 more states) at 
n = 9; and 39 more packings (and 70 more states) at n = 10. In a paper extending 
our work [30] . Hoy and O'Hern reported 52 packings at n = 9 and 279 packings at 
n = 10. However, the 279 packings that they reported multiply counted the three 25 
contact packings. When all instances of the same packing being counted more than 
once are removed, it yields 261 packings that were reported by Hoy and O'Hern 
in [30]. Please see supplementary information, Appendix B [8] for a list of exactly 
which packings are reported here but were not reported in [7j or in [30J. 

open square faces. We do not know whether or not at least 2 open square faces are a requisite 
of non-rigid packings that satisfy minimal rigidity constraints. The open square faces must be 
'connected' for the extra degree of freedom to exist - in the packings encountered thus far, this 
manifests itself by the existence of half-octahedra sharing at least 1 vertex. 

6.3. The Tree Nature of Packings. There is a distinct tree nature to packings. New seeds 
are the origin of a branch in the tree. Iterative packings continue the branch. All n particle 
iterative packings can be decomposed into combinations of less than n-particle packings, and this 
decomposition is often not unique. When the decomposition is not unique, the branches of the tree 
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Figure 13. New Seeds. 

All new seeds of n < 10 particles shown in both sphere and point/line represen- 
tation. There exists only 1 packing for each n of n < 5 particles, that can each 
be constructed iteratively from a dimer; thus there exist no new seeds for n < 5. 
n = 6 is the first instance of a new seed. The set of new seeds reported for n — 10 
is putative and thus represents a lower bound. New seeds with a * appearing to 
the right correspond to minima of the second moment out of all packings for that 
n. It can be seen here that, for the packings we have analyzed, when more than 1 
packing exists, the minimum of the second moment happens to correspond to a new 
seed. 



converge. Figures 16 and 17 show examples of tree convergence. Figure 16 shows the 2 possible 
decompositions of one of the 10 particle 25 contact packings. This packing can be formed either 
by adding 1 particle to the 9 particle non-rigid new seed, or by combining two polytetrahedra with 
a 6 particle octahedron. 

Figure 17 shows the tree structure of 2 < n < 8 packings. It can be seen that tree convergence 
occurs from n = 7 to n = 8, where multiple 7-particle packings produce the same 8-particle 
packing under the addition of another particle. 

6.4. Minima of the Second Moment. The second moment measures the deviation of particles 
from their collective center (centroid), and is given by 



(13) 
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(a) 




Figure 14. The 9 particle non-rigid new seed (grey) is shown in the top left-hand 
box. It is composed of two joined 5-point polytetrahedra (red) attached to two 
joined half octahedra (blue). The substructures are shown overlain - purple bonds 
and particles are shared by both substructures; whereas only red or blue bonds and 
particles belong to the polytetrahedral or octahedral structures alone, respectively. 
Two representative ways of forming this non-rigid structure are shown in (a) and 
(b). (a) Two 5-point polytetrahedra are joined by sharing one common point (on 
bottom). The two polytetrahedra are then fully attached to one another via the 
remaining 3 bonds first shown in dashed black lines, as potential bonds, and then 
in solid white lines, as actualized bonds. These bonds form the two connected half 
octahedra. (b) 2 particles (red) are attached to the concave side faces of the 2 joined 
half octahedra (blue). The 2 red particles form the two 5-point polytetrahedral 
substructures once they are attached to the joined octahedra. 

where rj is the x, y, z position of particle i; and c is the centroid, the average x, y, z position over 
all particles, given by 



and analogously given for c y and c z . 

The minimum of the second moment corresponds to the packing with the smallest M. The 2nd 
moment is listed within the list of packings in sectionJSl and in appendix B [8], and a '*' signifies 
the minimum of the 2nd moment for each n in figure [13] We confirm that the minima of the 2nd 
moment reported by Sloane et al. [48] are correct (they proved the 2nd moment minima for n < 4, 
but for n > 4 these were putative structures). For n < 10, each minimum of the second moment 
corresponds to a new seed, if a new seed exists. 

6.5. Ground States and the Maximum Number of Contacts. A fundamental question 
related to sphere packings is what is the maximum number of contacts that a packing of n spheres 
can have? Not only is this question of mathematical interest in its own right, but it is also 



(14) 




i=i 
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Figure 15. Non-rigid packings of 9 and 10 spheres, (a) The non-rigid n = 9 
packing, with non-rigid motion, corresponding to a twisting of the square faces, 
shown by black arrows, (b) Non-rigid n = 10 packings formed iteratively by adding 
one particle to the non-rigid n = 9 seed. The non-rigid motion of these structures is 
the same as in (a), (c) Non-rigid n = 10 seed. Non-rigid motion, corresponding to 
a twisting of the radially connected square faces, connected by the bottom particle, 
is shown by the black arrows. This twisting motion can be accomplished by, for 
example, twisting the top triangle. 



Figure 16. Tree Convergence in a 10-Particle Packing. 

An example of tree convergence in one of the 25 contact packings of 10 particles. 
This packing, shown in grey (top panel), can be decomposed into (bottom left 
panel) the 9 particle non-rigid new seed (red) plus one particle (blue) that rigidities 
the structure, or into (bottom right panel) an octahedron (blue) with 2 attached 
polytetrahedra (red). The red dashed line indicates an 'implicit contact,' a contact 
that automatically forms once the other contacts are in place (this corresponds to 
the 25th contact). 
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n = 2 






Figure 17. Tree Convergence for n < 8. 

An example of tree convergence for 2 < n < 8. Packings above which there is no 
arrow correspond to new seeds, and thus to the beginning of a new branch. The 
arrows then point to the n-particle packing or group of packings that form iteratively 
by adding one particle. It can be seen that even for the iterative case of adding one 
particle, there is tree convergence from n = 7 to n = 8 (shown by multiple arrows 
feeding into the same packing). 



of significant physical interest as such packings correspond to ground states. The number of 
packings that contain the maximum number of contacts in turn corresponds to the ground state 



degeneracy. Table 6.5 shows how the ground state degeneracy changes with n. Interestingly, this 
relation appears to be oscillatory. 

For n < 9, every packing has exactly 3n — 6 contacts. Thus, at each fixed n, n < 9, all packings 



have the same potential energy. Table 6J3 shows that, for n < 9, the ground state degeneracy 
increases exponentially. But at n = 10 this trend changes due to a small number of packings that 
can have 25 = 3n — 5 contacts (all other 10-particle packings have 3n — 6 contacts). There exist 
3 such packings, each containing octahedra (Fig. 18i-c). These three structures are the ground 
states at n = 10. 



NATALIE ARKUS, VINOTHAN N. MANOHARAN AND MICHAEL P. BRENNER 






(C) 






Figure 18. Maximum contact packings for 10 < n < 20. We are reasonably 
confident that the n = 10 packings shown here correspond to the maximum contact 
packings, but maximum contact packings of n > 10 are conjectured, (a)-(c) 10- 
particle packings with 3n — 5 = 25 contacts, (a) can be formed either by adding 
one particle (red) to one of the open square faces (blue) or as otherwise detailed in 
figure 16 (b) can be formed by adding one particle (red) to the concave 4 particle 
face created by 2 joined octahedra (blue), (c) is formed by connecting 2 octahedra 
by one edge, (d) 11 particle maximum contact packing (3n — 4 = 29 contacts). 
This can be formed either by adding one particle (red) to the concave 4 particle 
face created by the 2 joined octahedra of (b), or by adding one particle to the 
one remaining open square face of (a), (e) 12 particle maximum contact packing 
(3n — 3 = 33 contacts). This is formed by adding one particle (red) to the concave 
4 particle face created by 2 joined octahedra (blue). 13 particle maximum contact 
packings are constructed iteratively from this packing, (f) 14 particle maximum 
contact packing (3n — 2 = 40 contacts). This is formed by 3 radially connected 
octahedra (blue) and 2 particles (yellow) added to each of the concave 5 particle 
faces created by the 3 joined octahedra. (g) 15 and 16 particle maximum contact 
packings containing 3n — 1 and 3n contacts respectively. The 15-particle packing 
corresponds to the addition of only one of the red particles to the concave 4 particle 
face of (f), and the 16-particle packing includes both red particles, (h) 17 particle 
3n contact packing formed by 4 radially connected octahedra (blue) and 2 particles 
(yellow) connected to the concave 6 particle faces created by the 4 joined octahedra. 
At 17 particles, 3n contact packings correspond to this packing as well as packings 
constructed iteratively from (g). (i) 18,19,20 particle maximum contact packings 
corresponding to 3n + 1, 3n + 2, and 3n + 3 contacts respectively. Each of these 
packings is constructed by adding a particle (red) to one of the concave 4 particle 
faces created by the joined octahedra - the 18-particle packing is constructed by 
adding one such particle to (h), 19 by adding 2, and 20 by adding 3. 
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For n > 11, we conjecture as to what the maximum number of contacts are, and provide ex- 
amples of such structures. The maximum contact packings at n = 10 arise because it becomes 
possible to add 4 contacts to minimally-rigid 9-particle packings; whereas all other iterative pack- 
ings of n < 10 spheres are formed by the addition of 3m contacts to a minimally-rigid n — m sphere 
packing. All maximum contact packings found thus far correspond to iterative packings. We have 
not determined whether this is true for all n, but we conjecture that it is, because new seeds tend 
to contain more empty space, and thus less contacts. We have found three types of structures that 
allow for the addition of more than 3m contacts: (i) m octahedra, where each pair of octahedra 
share one edge (as in Fig. [l8b), (ii) an open square face created by half an octahedra (as shown 



in blue in figure 18 1), and (iii) the concave m point face created by octahedra sharing 3 edges (as 
shown in blue in figure 
concave faces in figure 



18i,e). 4 point concave faces are shown, for example, in figure 18 3; 5 point 
18f,g; and 6 point concave faces in figure |i~8~p,i. To each m point concave 
face, it is possible to add one particle with m contacts; this is evidenced, for example, by the red 
particle with 4 contacts in figure |l8fe and the yellow particles with either 5 or 6 contacts in figure 
18~j>i. All structures leading to maximum contact packings that we have found thus far are related 
to octahedra, and we conjecture that this will be the case for all n. (Interestingly all non-rigid 
packings found thus far are also related to octahedra, in that they contain the open square faces 
created by half-octahedral structures; they all contain half-octahedra sharing at least 1 point). 

There are many fewer ways of adding greater than 3m contacts to an n — m sphere packing 
than there are of adding 3m contacts, thus leading to a relatively small number of ground state 
packings when a new maximum number of contacts, as a function of n, is reached. Thus, each time 
a packing with a greater maximum number of contacts, as a function of n, is possible, we expect 
the ground state degeneracy to either decrease or to remain small. When the functional form 
for the maximum number of contacts remains constant, we expect the ground state degeneracy 
to grow rapidly, due in large part to the iterative growth of adding a particle with 3 contacts 



to packings. Table 6.5 for example, shows that for n < 9 the ground state degeneracy increases 
exponentially because all packings have 3n — 6 contacts. But at n — 10 this trend changes because 
packings with more than 3n — 6 contacts become possible. At n — 13 and at n = 17, rapid 
growth in the ground state degeneracy resumes as the functional form for the maximum number 
of contacts remains constant, but for all other n, the ground state degeneracy either decreases or 
remains constant, because a new maximum number of contacts, as a function of n, is possible. 

At low temperatures, we expect the experimentally observable packings to be dominated by the 
maximum contact packings. Thus, for n = 10, and for any n in which either (i) the maximum 
number of contacts increases as a function of n or (ii) the maximum number of contacts does not 
remain constant for too long (< 3 particles), we expect that there are only a small number of 
packings that will be observable at low temperature. This trend can be seen in the conjectured 



ground state degeneracy for 9 < n < 20 in Table 6.5, and we conjecture that a similar low 



amplitude oscillatory trend might also continue for n > 20. 

6.6. Lattice Structure. The non-rigid new seeds at n = 9 and n = 10, as well as the maximum 
contact packings of n < 13 are all sub units of the hexagonally close-packed (HCP) lattice, being 
combinations of face-sharing tetrahedra and octahedra. Additionally, the structure shown in figure 
is a subunit of either the HCP or face-centered cubic (FCC) lattice. The non-rigid packings 



are entropically favored, and we thus expect these to form with higher probability at higher 
temperatures; while the maximum contact packings are energetically favored, corresponding to 
the structures that will form with higher probability as the temperature, T — > 0. For 14 < n < 20, 



the maximum contact packings are not HCP subunits (Fig. 18 f-i 
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Table 5. Number of Contacts vs. Number of Particles. 

The first column corresponds to the number of particles, n, the 2nd column to the 
number of contacts that the ground state packing(s) have, and the 3rd column 
to the number of packings in the ground state. Note that for n > 11 these are 
putative results and thus the number of packings with 3n — m contacts represents 
a lower bound, and the ground state number of contacts is conjectured (we have 
not encountered packings with more contacts, but we have not proved that they 
don't exist). 



n 


L-ontacts 


Ground State Packings 


2 


3n — 6 


1 


3 


3n — 6 


1 


4 


3n — 6 


1 


5 


3n — 6 


1 


6 


on — 6 


I 


7 


3n — 6 


5 


8 


3n — 6 


13 


9 


3n — 6 


52 


10 


3n — 5 


3 


11 


3n-4 


1 


12 


3n — 3 


1 


13 


3n — 3 


7 


14 


3n-2 


1 


15 


3n - 1 


1 


16 


3n 


1 


17 


3n 


8 


18 


3n + 1 


1 


19 


3n + 2 


1 


20 


3n + 3 


1 



Frank predicted [20] that icosahedral short-range order would be a hallmark of liquid structure, 
and experimental studies have shown local cluster-like order in bulk atomic liquids and glasses 
[HI ST] . Results from a recent study suggest that structural arrest in condensed phases may be 
related to geometrical constraints at the scale of a few particles |45j . The propensity for icosahedra 
[28| 116] in longer-range systems is absent in ours. We have proven that the icosahedron is not the 
ground state at n = 12, nor is an icosahedron with a central sphere the ground state at n = 13. 
A 12-sphere icosahedron has only 3n — 6 = 30 contacts, and in a 13-sphere icosahedron the outer 
spheres would not be close enough to interact with each other. 

It is possible, and perhaps even likely, that the lattice structures corresponding to ground state 
packings will be periodic with n. For example, although the ground states for 14 < n < 20 are 
not commensurate with HCP, the ground states for a finite range of higher n may be, and may 
then subsequently return to the lattice structure commensurate with 14 < n < 20. Detailing 
the structures of ground state packings for all n, and geometrical patterns contained therein, is a 
subject of future work. Furthermore, the appearance of crystalline order, such as HCP, at very 
low n may influence nucleation. 
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7. Extensions and Conjectures 

7.1. The major roadblock for reaching higher n. The main roadblock to the analytical enu- 
meration of sphere packings at higher n in the current work is deriving one analytical geometrical 
rule that can solve for all new seeds. In the next section, we outline a numerical method, based 
on the triangular bipyramid rule, that is capable of finding all solutions of A — > V, and which 
can thus solve for all new seeds. However, the implementation of either this numerical scheme 
or the derivation of an analytical rule would only allow us to enumerate packings of up to about 
n = 14 spheres. This is because the real limitation of the current method arises from the enu- 
meration of minimally-rigid, non-isomorphic adjacency matrices. For n < 10, the enumeration 
of such adjacency matrices using nauty [12] takes on the order of seconds. For n = 10, 11, this 
enumeration takes minutes. For n = 12, enumerating all minimally-rigid, non-isomorphic .4,'s 
takes approximately 2 hours. Extrapolating, we expect the enumeration at n — 13, 14 to take 
on the order of 2 days and 2 weeks, respectively. Thus, around n = 14, we begin to reach the 
computational limitations of this method, which is due to the enumeration of A's. 

Only a very small fraction of adjacency matrices correspond to sphere packings; for example, at 
n = 10, out of the 750,226 .A's, only 262 correspond to sphere packings. Thus, the enumeration of 
all A's really is a brute force and wasteful step. Further advances in enumerating sphere packings 



will require overcoming this roadblock. In section 7.3, we propose one method that might be able 
to overcome this limitation. 

7.2. Applying the Triangular Bipyramid Rule to New Seeds. The Triangular Bipyramid 
rule solves for iterative structures but does not work for noniterative ones, which we also showed 
increase rapidly starting at n — 10. Here, we discuss how the triangular bipyramid rule might 
also be applied to new seeds. In this case, the equations for the unknown interparticle distances, 
rjj, are implicit and thus must be solved numerically. 

In a new seed, any triangular bipyramid that contains an unknown distance will contain more 
than one unknown distance. Thus, can not be determined one at a time, as with iterative 
packings; instead we must construct a system of equations to solve for the unknown distances; 



Let us consider a general triangular bipyramid, like that shown in figure 19, but here all 10 



distances are potentially unknown (Fig. 19). For m unknown r^, we construct m triangular 
bipyramids to solve for the (see Fig. 20, for example). Each triangular bipyramid yields one 
equation for an unknown distance, thus yielding m equations with m unknowns in total, making 
the system well defined. 



We assign a label, di, to each of the 10 distances within the triangular bipyramid (Fig. 19). 
Explicit formulas can always be obtained for d 2 and d%. Thus, over all m triangular bipyramids, 
we place each unknown distance in location d 2 or d 3 at least once. New seeds are inherently 
global structures, thus the m triangular bipyramids should contain all n points amongst them in 
order to ensure that the solution space is sufficiently constrained. Also, to avoid redundancy, each 
triangular bipyramid must contain a unique combination of 5 particles. 

The equations for the fy are derived from the triangular bipyramid property that the dihedral 
angle A\ must equal either the sum or the difference of the dihedral angles A 2 and A% (as was 



detailed in section 4 and in Fig. 19). This equation is cumbersome, and it as well as its derivation 



can be found in Appendix A of the supplemental information |8J. As always with the triangular 



19 As can be seen in section [4] or in Appendix A of the supplemental information [8], the equation for an unknown 
distance is quadratic. Thus, this will be a quadratic system of equations. 
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Figure 19. General Triangular Bipyramid. A triangular bipyramid where 
all distances, labeled d\ — d w , are potentially unknown. The particles, i,j,k,p,q 
correspond to the same labeling as in figure [T9j This triangular bipyramid is used 
to derive the general equation for an unknown relative distance, (appendix A 

0)- 



bipyramid rule, due to the 2 possibilities of A\ equaling either the sum or the difference of A 2 and 
A 3 , each unknown distance has 2 possible solutions. 

For each set of that is to be solved, construct initial guesses between the bounds of the 
triangle inequality and no-overlap constraint, and iterate the initial guess with a step size less 
than or equal to the minimum difference between different solutions (for rigid structures). There 
will always exist unknown < 2R because each particle has at least 3 contacts^ These are 
the Tij for which there exists a k satisfying Aik = Ajk = 1, A\j = 0. Thus, first solve the set of 
R < r ij < 2i?. If unknown remain, then solve the set of that now have known triangle 
inequality bounds, due to the previously solved set of r^, repeat until all have been solved^ 1 ] 

7.3. The Bond Breakage Conjecture. A packing of n spheres can be formed by (i) taking 
an n — m sphere packing, breaking a contact (or bond), adding m new spheres, and forming the 
appropriate contacts to complete the packing, or by (ii) breaking one bond of an n sphere packing 
and reforming another. From this property, we propose the following theory. 

Bond Breakage Conjecture: All packings of n spheres can be obtained by breaking one bond 
and reforming another in every possible way. For any packing, there exists an m step path, of 
breaking one contact and reforming another, that will form that packing out of another packing 
with 3n — 6 contacts. Each of the m steps will correspond to an n-particle packing. The end points 

20 This also holds true when each particle has at least 2 contacts. 

21 We have tested this method on the new seed ,4's for n < 8, and have shown that it works; however, we have 
not implemented it for up to n = 14. 
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Figure 20. A New Seed and the General Triangular Bipyramid. 

Here we show an example of how to apply the general triangular bipyramid to a new 
seed in order to determine its unknown distances. This new seed corresponds to the 
6 particle octahedron. It has 3 unknown distances, r\ 2 , r^, and (dashed black 
lines). The first triangular bipyramid constructed consists of particles 1,2,3,4,5, 
and has unknown distances d\ and d 2 , corresponding to r 34 and r 12 , respectively. 
Note that (d 3 to d w = R). The second triangular bipyramid consists of particles 
1,2,3,4,6 and has unknown distances d' 1 ,d' 2 corresponding to r 12 ,r 34 , respectively. 
The third triangular bipyramid consists of particles 1,2,4,5,6 and has unknown 
distances d", d 2 corresponding to r\ 2 , r^, respectively. Note that the first two trian- 
gular bipyramids comprise 2 equations and 2 unknowns, and thus are alone sufficient 
to solve r\ 2 and r^. Once these 2 distances are known, applying the third triangular 
bipyramid involves only 1 unknown distance, r 56 , and thus follows as applying the 
triangular bipyramid to an iterative packing. 

of the path (i.e. which 3n — 6 packing one begins with and which packing one ends up with) deter- 
mine what value m takes. For every packing, there exists at least 1 other packing for which m = 1. 

This suggests an alternative method for enumerating all packings of n spheres: construct just 
one 3n — 6 contact packing of n particles (this can easily be done, simply construct an n particle 
polytetrahedron for example), and then break and reform bonds in all possible ways. For each 
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packing, it is important to explore every combination of breaking and reforming a bond - i.e. to 
go down all paths, and not just one path. 

We have confirmed this conjecture up to as high as we have enumerated packings (n = 10) 
using the following algorithm: 
For every A that corresponds to a packing 
1. For each 1 that appears in the A: 

(a) Swap the 1 with an existing 0. Do this in every possible non-isomorphic way. (This 



is the mathematical analogue of physically breaking an existing bond and reforming 
a different bond that was not present in that packing). 

For each new A that is generated by swapping a 1 and a (these are 1 bond away 
A's, i.e. where m = 1), 

(i) Test for an isomorphism with the A's of all other packings (i.e. all .A's other 
than the one being examined). 

(ii) If an isomorphism is found, stop the examination of this A, as it has been shown 
that there exists a 1 bond bath between the packing being examined and another 



Thus implementing this algorithm computationally, we have proven that, for each n < 10, 
every packing is a 1 bond distance away from at least 1 other packing of the same n. We have not 
proven this for n > 10, as we have not enumerated all packings of n > 10, but we suspect that 
this conjecture holds for all n. 

Mapping out all possible 1 bond distances might be able to shed insight into the kinetic pathways 
of packings. 

To implement the bond breakage conjecture into an improved method for enumerating sphere 
packings, the bonds must be broken and reformed intelligently, such that unphysical conformations 
are not explored. If all physical and unphysical conformations were to be explored, then we would 
simply return to the same computational problem we had with enumerating all non-isomorphic 
.A's. One should be able to break a bond and reform it only as is physically possible by calculating 
the 1 degree of freedom motion that is left over from the one broken bond. Furthermore, one can 
take advantage of symmetry to, a priori, not explore redundant (i.e. isomorphic) pathways of 
breaking and reforming bonds. 

7.4. Extensions to Other Dimensions. The method presented here is, in large part, not di- 
mension specific. The only step which is dimension specific is the set of geometrical rules used 
to solve A for T>. However, at least for the triangular bipyramid rule (which solves for most of 
the packings), the geometrical rules can easily be modified to account for a different number of 
dimensions, d. Once this is done, the same method can be used to solve for sphere packings of 
d ^ 3 dimensions. 

7.5. 'Lower Dimensional Packings'. While a packing of n spheres depends on the dimension- 
ality of its embedding space, there exists a cutoff number, m, for which the packing of n spheres 
remains constant for all d + i dimensions, i > 1. For this m, the n spheres have accessed all 
dimensions possible to them, and so the embedding space becomes irrelevant. For example, for all 
d, the unique packing of 1 and 2 spheres is the singlet and doublet, respectively. For 3 spheres, the 
unique packing in 1 dimension is a linear connected chain of 3 spheres; whereas in 2 dimensions it 




If one is interested in examining all 1 bond paths that exist between all packings, then this same algorithm 
can be executed without this termination step to yield all possible 1 bond paths. 
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is the equilateral triangle. For d > 2 however, the unique packing of 3 spheres remains the same, 
it is always the triangle. For 4 spheres, the unique packing in 2 dimensions is the di-triangle; 
whereas, in 3-dimensions it is the tetrahedron. It is generally true that packings of d + 1 particles 
remain the same for d or more dimensions^ 3 ! 

7.6. Patterns in Adjacencies and Distances. Does there exist a signature pattern in the ,4's 
orD's that signifies a packing? In other words, is there a pattern in the distribution of adjacencies 
(i.e. number of contacts per particle and connections therein) and/or the distribution of distances 
that corresponds to packings? If such a pattern does exist, it could illuminate a more general 
method for solving for packings. It also might shed light on the spectrum of allowed solutions for 
the system of quadratic equations corresponding to the adjacency matrices - detailing which do 
and do not have real valued solutions in IR 3 satisfying T> > R. 

7.7. Related Mathematical Problems. 

7.7.1. Erdds Unit Distance Problem. The Erdos unit distance problem (a.k.a. Erdos Repeated 
Distance Problem)^] was posed in 1946 by the Hungarian mathematician, Paul Erdos. It asks what 
the maximum number of unit (or repeated) distances that can connect n points in d dimensions 
is [121 E]. This problem is still unsolved. Even in 2 and 3 dimensions, only upper and lower 
bounds are known [TU [22] • The solution to this problem in 3 dimensions, where the unit distance 
is also the minimum distance, would answer what the maximum number of contacts in any sphere 
packing is; thus giving the number of contacts corresponding to the ground state packing(s). 

7.7.2. 3-Dimensional Rigidity. In solving adjacency matrices for both rigid and non-rigid pack- 
ings that satisfy minimal rigidity constraints in 3-dimensions, this problem is directly related to 
determining whether a graph is rigid in 3-dimensions. Much work has been done in this field 
[33] [13] 135] [12] [32], [22], as well as in other dimensions [40]. The existing work on rigidity may help 
to further inform sphere packings, and the work presented here may in addition be applicable to 
rigidity theory. In particular, it may allow for the development of a simple method for reading off 
whether a 3-dimensional graph is rigid or not. By the method presented here, a graph is deter- 
mined to be non-rigid if there exists a continuum of solutions to A. However, if we can determine 
a signature pattern that corresponds to all non-rigid (but minimally-rigid) A's, this would allow 
for a very simple determination of whether a 3D graph is rigid. 

7.7.3. Solutions to Systems of Polynomial Equations. The method presented here is inherently 
solving a system of quadratic equations. Thus presenting an alternative analytical solution to this 
class of problems. Current standards in the field for analytically solving systems of polynomial 
equations include Grobner bases [TT]; however, these are time-consuming and thus do not scale 
efficiently with the number of equations. The method presented here solves a certain class of 
polynomial equations efficiently for a relatively large number of equations. Is it possible to extend 
this method in order to more efficiently solve large systems of polynomial equations? 

23 Wc originally posed this as a question, but one of the reviewers of this paper pointed out that this was true 
due to the fact that any d + 1 points are always contained in an affine subspace of at most d dimensions (without 
loss of generality, one point can be taken to be the origin). Thus, the d+ 1 points can only access d dimensions, 
and considering a higher than d dimensional space will not change the structures that they can form. 

24 Without loss of generality, a repeated distance can be called a unit distance, because one can always uniformly 
rescale all distances such that the repeated distance is the unit distance. Put another way, 'a unit' can be given 
any value - here, the unit is simply given the value of the repeated distance. 
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7.7.4. Euclidean Distance Matrix Completion Problems. Given a symmetric matrix, M, where 
only certain elements are specified, the Euclidean distance matrix completion problem is to find the 
unspecified elements of M that make M a Euclidean distance matrix. Euclidean distance matrix 
and positive semidefinite matrix completion problems are closely linked [361 E3 EEJ EH EU Hj. In 
solving adjacency matrices for distance matrices, the method presented here is directly related 
and potentially directly applicable to the euclidean distance matrix completion problem and, by 
extension, to the positive semidefinite matrix completion problem. 

8. Concluding Remarks 

In this work, we present an analytical method for deriving all packings of n spheres. We carry 
out this derivation for n < 10; where the set of n = 10 new seeds is preliminary, and all iterative 
packings of n = 10 spheres and all packings of n < 9 spheres are potentially complete, save the 
numerical round-off error present from implementing this analytical method computationally. 

We consider the derivation of these sphere packings to be the first step in directing the self- 
assembly of spherical colloidal particles; where we have divided this problem up into 2 parts: (i) 
understanding what the system of colloids can self-assemble, (ii) deriving a mechanism to control 
that self-assembly. The derivation of all packings of n spheres gives us everything that a system 
of n colloidal particles can self-assemble, thus taking care of this first step. Future work will 
detail the second step, which is the derivation of a mechanism that directs the self-assembly of 
the packings such that only one packing forms. 

Beyond the problem of self-assembly, the results reported here are interesting in their own right. 
We find many interesting properties from the sphere packings enumerated up to n = 10, as well 
as from the conjectured maximum contact packings of 11 < n < 20. Furthermore, the results are 
directly related to the physics of colloidal clusters, and may have applications to glassy systems 
and the nucleation of crystals. They are also directly related to unsolved problems in mathemat- 
ics, such as the Erdos unit distance problem. 

We thank John Lee for consultations in coding, and for writing code that (i) converted the pack- 
ing output to latex format for appropriate display in appendix B, and (ii) interfaced the packing 
output with the rotational constant calculator website and with the mayavi2 graphing package so 
that the 10 particle point groups, symmetry numbers, and figures did not have to be generated 
manually. We also thank David Roach, Noam Elkies, and Marcus Roper for helpful discussions. 
We acknowledge support from the MRSEC program of the National Science Foundation under 
award number DMR-0820484, the NSF Division of Mathematical Sciences under grant number 
DMS-0907985, and DARPA under contract BAA 07-21. 
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APPENDIX A: INDIVIDUAL GEOMETRICAL RULES 



In this appendix, we derive the individual geometrical rules that were used to solve A for T>. 
1. Rules Sufficient to Solve All 6 Particle Packings 

Rule 1: 

The first rule is the most simple; it states that any 2 intersection circles can not intersect in 
> 2 points (fig 2). 




(a) 2 intersections of inter- (b) 1 intersection of 2 intersec- (c) intersections of intersection 
section circles tion circles circles 



Figure 1. Intersecting Intersection Circles 

This rule is identified by the following adjacency matrix pattern: 

Any 2 of the set {Ajk, Aj P , Ak P } equal 1, and 3 > 2 z's for which Aij = Aik = Ai P = 1. 

Rule 2: 

The 2nd rule calculates the distance between the 2 points of intersection of 3 mutually intersect- 
ing intersection circles associated with a triplet (fig [2j). If a particle lies at one of the intersection 
points, this forms a 4 particle packing (the tetrahedron), for which all distances = R. If 2 particles 
lie at each of the intersection points, this forms a 5 particle packing - the distance between these 
2 points is the only distance > R. 

The non-unit distance (where the 'unit distance' is = R) of the 5 particle packing is calculated 
as follows (fig 4): A trimer forms an equilateral triangle with sides of length R. The intersection 
circles have a radius of (\/ r 3/2)R, and the points, corresponding to the particles' centers, lie on 
the circumference of these intersection circles. The centers of the intersection circles lie at the 
midpoint of each of the triangle sides (of length R). Thus, we can draw a line segment of length 
(v3/2)i2 emanating from the points of the equilateral triangle to the midpoints of the opposite 
line segment. (See figure [3]). The intersection of these line segments corresponds to the in-plane 
location of where the intersection circles intersect - i.e. projecting the intersection points to the 
face of the equilateral triangle. For example, in the pentamer in figure [3|d (i.e. the 5 particle 
packing), the center triangle corresponds to the in plane triangle we are considering, and the top 
and bottom points correspond to the 2 intersection points of the 3 intersection circles. As all 
of these line segments = R (as they represent particles touching each other), all of the triangles 
formed in this pentamer are equivalent equilateral triangles oriented differently in 3-dimensional 
space. Let us consider the 3rd triangle in figure [3^,; using right triangle math, we can calculate 
what a is: 
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4 particle packing: 




5 particle packing: 




Figure 2. The Intersections of 3 Intersection Circles 



tan 



7T 



6 
1 

71 



IR 



t;R 



2^3 



R 



Because all of the triangles in the pentamer are equivalent, we can use a to construct the triangle 
in figure where the (y/3/2)R edge corresponds to the radius of one of the in-plane intersection 
circles. Thus, the (\/3/2)R edge connects one of the in-plane triangle edges to the point of 
intersection of the 3 in-plane intersection circles. The l/2h edge corresponds to the line segment 
connecting the point of intersection of the 3 circles of intersection to the in-plane projection of this 
intersection point; and therefore represents 1/2 the length between the 2 intersection points (as 
the in-plane trimer lies at the midpoint of the 2 intersection points). This forms a right triangle, 
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V3/2 R - a 





1/2 h 



1/(2 V3) R 



Figure 3. Calculating the Distance Between the 2 Points of Intersection of the 3 
Intersection Circles of a Trimer (or Triplet). 



and once again we can use right triangle math to solve for h: 



1 



2V3 



R) + 



-h 



x/3 



R 



(15) 



h 



2 < R 



Thus, the solution to an adjacency matrix corresponding to a 5 particle packing, written in 
terms of a distance matrix, is 

/Oil 1 1 \ 
10 11 1 
110 1 1 



/0 1 1 1 l\ 
10 111 
110 11 
1110 

\1 1 1 0/ 



1 1 1 
1 1 1 







/ 



This rule can be identified by the following adjacency matrix pattern: 



A 



1, and 3 2 p for which A P i = A P j = A p k = 1. This rule solves for the distance 



between those 2 p, for ex. if p = I, m, then V 



lm 



TO 

2JjR. 
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Figure 4. 4 Touching Points on an Intersection Circle 



Rule 3: 

This follows from rule 2, namely as we have calculated that the distance between the 2 inter- 
section points of 3 mutually intersecting intersection circles is 2^j2/3R, then that distance can 
not = R. Any adjacency matrix implying that distance = R is unphysical. 



This rule can be identified by the following adjacency matrix pattern: 



A%j A, 



if. 



A, 



1, 3 2 p for which Am = A 



X P3 



and the A element associated with those 2 p's is 1 



A p k = 1 (he. the same pattern as in rule 2), 
[i.e. if p = q, z then A qz = 1). 



Rule 4: 

This rule calculates the distance between the two end particles of a total of 4 particles that 
lie on a circle of intersection, where particle i + 1 touches particle i (1 > i < 4). (See figure [4]). 
From the figure, we can see that this scenario forms 3 equivalent isoscoles triangles, which we will 
be able to use to calculate the distance, r, between the 2 end particles (see figure [5]). We can 
calculate 9 from this triangle, either via the law of cosines 

( (fR) 2 + (fi?) 2 - R 2 \ 



cos 



V 



2^R^R 



J 



cos 



-1 



or via right triangle math: 



IR 
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Figure 5. Solving the Long Distance Between 4 Touching Points on an Intersection Circle 

and both of these expressions for 9 are equivalent (both « 1.23096). 

Using 9, we can calculate the desired distance, r. The new angle, 0, depicted in figure [5] is given 
by 2vr - 39: 



and thus r is given by 



and as 




r is given by 

(16) r = h -R 

Therefore, the distance between particle 1 and 4, where particle i + 1 touches particle i for i — 1 
to i — 4, and particles 1 to 4 all touch the same 2 particles is always 5/3R. 
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1 
1 








r 
















Figure 6. 4 Intersection Circles Intersecting at 2 Points 

The red points in the first figure are the points of intersection of the 4 intersecting circles. The 
dashed parts of the circle indicate that part of the circle goes into the page, while the solid lines 
come out of the page. The second figure, drawn with one intersection circle not shown, is added 
for clarity, so that the intersection points of the intersection circles can become more apparent. 



This rule can be identified by the following adjacency matrix pattern: 

Aij = 1, and 3 4 for which Aik = Ajk = 1, and 3 A pq = 1 amongst the 4 k (i.e. if k — p, q, I, m, 
then A pq = A q i = Ai m = 1), and the distance T> pm = 

Rule 5: 

Following from rule 4, we know that any adjacency matrix implying that the above distance 
= R is unphysical. 



This rule can be identified by the following adjacency matrix pattern: 
The same pattern as in rule 4, where A pm = 1 also. 



Rule 6: 

This rule derives the distance between the 2 intersection points of the mutual intersection of 4 
intersection circles associated with 4 points (fig [6]). As mentioned before, m intersection circles 
can intersect at 1, 2, or points. As we can see from figure [6j when 4 intersection circles intersect 
at 2 points, they must be associated with a 2-dimensional square (in fact, we will see later that 
this is generally true - if it is possible for m intersection circles to intersect at 2 points, they will 
only intersect at 2 points when associated with a regular m-gon lying in a 2 dimensional plane. 
And 2<m<5, ifm>6 the m associated intersection circles can intersect in at most 1 point (see 
rule 13). If 4 intersection circles are not associated with a square, then they can only intersect at 
1 or at points. 

We can calculate the distance between the 2 points of intersection of 4 mutually intersecting 
circles by constructing an isosceles triangle that connects one of the intersection points to the 
center of 2 of the intersection circles - such a triangle is shown in figure [7j We know the lengths 
of all of the sides of this triangle, as the 2 equivalent lengths are the radius of the intersection 
circle, (v / 3/2)i?, and the base length is the distance between 2 touching particles (or equivalently 
the radius of the neighbor sphere), R. The height of this triangle corresponds to half the distance 
between the 2 points of intersection, thus we can use right triangle math to calculate the desired 
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V3/2 R 



1/2 h 




1/2 R 



Figure 7. Solving for the Distance Between 2 Intersection Points of 4 Mutually 
Intersecting Intersection Circles Associated with 4 Points. 



distance: 



4 4 



V3 



1? 



4 4 



Z? 



(17) ==> [fc = y^i? 

Thus, the distance between 2 points of intersection of 4 intersection circles must be V%R, if any 
adjacency matrix violates this property, by setting this distance = R, then it is physically impos- 
sible. 



This rule can be identified by the following adjacency matrix pattern: 

3 a closed 4 ring, A%j = Ajk = Ak p = A p % = 1, where the cross particles don't touch, Aik = Aj P = 
0, and 2 particles touch all points in the 4 ring, A m j = A m k = A mp = A m i = 1 for two m's. This 
rule then solves for the distance between those 2 m's, i.e if m = z, s then V zs = y/2R. 

The following adjacency matrix contains this pattern: 

/0 1 1 1 1\ 

1111 

110 11 

110 11 

11110 
\1 1 1 1 0/ 



Rule 7: 

As just mentioned, following from rule 6, any adjacency matrix that implies the distance be- 
tween the 2 intersection points of 4 mutually intersecting intersection circles = R is unphysical. 
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s 




Figure 8. 5 Points on an Intersection Circle. (The '?' distance = the unknown 
distance r.) 

This rule can be identified by the following adjacency matrix pattern: 

The same pattern as in rule 6, and the adjacency matrix element between those 2 m's is 1 (i.e. if 
m = q,l, then A q i = 1 also). 

2. Rules Sufficient for 7 Particle Packings 

The following geometrical rules, in addition to the rules used to solve the 6 particle packings, 
are sufficient to solve for all 7 particle packings. 



We have calculated the long distance between 4 and between 3 touching points on an intersection 
circle (rules 2 and 4); for 7 particle packings, we must also know the long distance between 5 
touching points on an intersection circle (i.e. the distance between points 1 and 5, where point 
i + 1 touches point i for 1 < i < 5) - see figure [8] 

We make this calculation as we did for the one with 4 touching points on an intersection circle, 
but here the relevant angle of the triangle that will allow us to calculate the distance r is 



Rule 8: 




2n — 8 sin 
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and thus r is given by 



sm 



sm 7T 



4 sin 



-i 



(18) 



-6 



1 

Ay/2 
~9~ 



R 



v/3 



R 



R 



R 



This rule can be identified by the following adjacency matrix pattern: 

Aij = 1, and 3 5 A; for which A^ = Ajk = 1, and A A pq = 1 amongst the 5 k (i.e. if k — p, q, I, m, z, 
then A pq = A q i = Aim = A mz = 1). This rule solves V pz = ^-R. 

Rule 9: 

Following from rule 8, we can see that 5 particles can not touch the same 2 particles and all 
touch each other, as the distance between the 4th and 5th points lying on an intersection circle is 
= ^-R ~ 1.08866-R, and for 5 particles to touch the same 2 particles, this distance would have 
to equal exactly R. Thus any A implying this distance = R is unphysical. 



This rule can be identified by the following adjacency matrix pattern: 

The same pattern as in rule 8, but 5 A pq = 1 amongst the 5k- i.e. A zp = 1 also. 

Rule 10: 

Also following from rule 8, we can see that no more than 5 points can fit on an intersection 
circle (i.e. no more than 5 points can mutually touch a dimer), as only slightly more than 5 points 
can fit on an intersection circle, and 6 points can not fit. 

We can calculate exactly how many points a distance R apart can fit on an intersection circle 
by dividing the circumference of the entire intersection circle by the arc length swept out by one 
of the isosceles triangles created by 2 particles lying on the intersection circle. (See figure [8]). The 
formula for the arc length is 

S = r8 

where r is the radius of the circle, 9 the angle between 2 radial line segments, and S is the arc 
length. Thus the number of points a distance R apart that can fit on an intersection circle is 

2tt^R tt 

^008-1 (|) " COS-1 (!) 

(19) ss 5.1043 
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Figure 9. Calculating the Distance Between 2 Intersection Points of 5 Mutually 
Intersecting Intersection Circles Associated with 5 Points. 



(As an aside, cos _1 (l/3) = 2 sin _1 (l/\/3) = 2 tan _1 (v^2/2), and we could have expressed the angle 
by any one of these). 

Thus any A implying that > 5 points lie on an intersection circle is unphysical. 



This rule can be identified by the following adjacency matrix pattern: 
Aij = 1, and 3 > 5 k for which Aik = Ajk = 1. 



Rule 11: 

As we did previously with the mutual intersection of 4 intersection circles (rule 6), we will 
now calculate the distance between the 2 intersection points of 5 mutually intersecting intersec- 
tion circles. 5 circles of intersection can only intersect at 2 points if they are associated with a 
2-dimensional pentagon (see figure [9]). To calculate the distance between these 2 points of inter- 
section, we can construct the triangle that connects one of the points of intersection to the centers 
of one of the intersection circles (shown in the figure). The base of this constructed triangle, a, 
is the apothem of a regular pentagon. We can calculate the apothem of the regular pentagon via 
the right triangle depicted in figure 10 (the formula for the apothem is also well known, so one 
could alternatively simply look it up). The interior angle of the pentagon, <p (fig 10), is given by 

(n - 2)vr 3 

(f) = = -7T 

n 5 

and the angle associated with the right triangle used to calculate the apothem is (f>/2. Therefore, 
the apothem, calculated via right triangle math, is given by 

, 3 

tan — 7i 



10 J \R 



1 + R 
2^10-2^ 
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Figure 10. The Apothem and Interior Angle of the Pentagon. 



Now that we know a, we can calculate the distance between the 2 points of intersection, h, via 
right triangle math: 



(20) 



a 2 + \-h 



] + VE r) +\« 

,2^10-2^ / 4 
ll + 2V5 + 5 fl2+ l ft2 
4 10-2^ 4 



'VI 

2 
l R 



R 




This rule can be identified by the following adjacency matrix pattern: 

3 a closed 5 ring, Aij = Ajk = Ak p = A pz = A Z i = 1, where the cross particles don't 
touch, Aik = Aj P = Akz = A P i = Aj z = 0, and 2 particles touch all points in the 5 ring, 
A m j = A m k = A mp = A mz = A m i = 1 f° r two m's. This rule identifies the distance between the 

2 m's (i.e. \im — l,y this rule gives T>i y = J2 — -^R)- 
Rule 12: 

Related to rule 11, the distance between the vertices of the pentagon will be the distance between 
the non-touching pairs amongst the 5 particles that have intersection circles mutually intersecting 
at 2 points (X^, T>k z , T> pi , T>j Z from rule 11). This distance can be calculated from the triangle 



in figure 11 We can solve for the distance, d, by using right triangle math or by using the law of 
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Figure 11. A Triangle Within the Pentagon 



cosines, here we'll use the law of cosines 

d = 



R 2 + R 2 - 2R 2 cos [ -7T 



'R 2 + R 2 - 2R 2 (-^ (-1 + Vtj^j 



(21) 



d 




This rule can be identified by the following adjacency matrix pattern: 

The same pattern as in rule 11, but this rule identifies the distances between the cross particles; 



V 



Rule 13: 



V kz = V. 



pi 



V 3Z 



R 



6 or more mutually intersecting intersection circles can only intersect at one point. A regular 
m-gon is the limiting case where m intersection circles will intersect at 2 points, if they can. For 
3, 4, and 5 particles, this is possible. For > 6 particles however, this is not. The longest distance 
between the non-touching particles of a regular hexagon, whose sides are length R, is 2R. (See 
figure 12). The interior angle of a regular n-gon is 

(n - 2)tt 



n 



Thus, 6 here is 



-7T 



The angles associated with the equivalent sides of the isosceles triangles formed by the n-gon is 
((n — 2)/2n)7r, which implies that 6 is given by 

(n - 2)tt 







7T — 



n 
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Figure 12. The Hexagon 



see fig 12 ) and we can see that h is indeed equal to 2R, from right triangle math (we could also 



have done this, as usual, from the law of cosines): 



sin —9 



l„\ hR 



2_ 

2) \h 
1 \ R 

Sin I -7T 



6 J h 

1 R 

2 ~ ~h 

(22) => h = 2R 

This distance is equal to twice the distance of 2 touching particles (i.e. to the diameter of the 
neighbor sphere), and thus we can see that one particle fits right at the center of the regular 
hexagon, showing us that 6 mutually intersecting intersection circles can only intersect at one 
point (as only 1 point can touch all 6 points of the hexagon). 

Put another way, we can calculate that the height of one of these isoscoles triangles is = 
(y3/2)R, thus implying that the distance from the center of an intersection circle associated with 
one of the hexagon edges is a distance a/3-R away from the center of an intersection circle associated 
with the opposite edge. Therefore, the distance between the center of the intersection circle and 
the intersection point of the intersection circles is (v / 3/2)i?, and thus the 'height between the 2 
intersection' points must be 0. On the other hand, we can see that the distance between the 
longest non-touching vertices of a regular septagon is > 2R (not shown here), and thus 7 in-plane 
intersection circles can not all intersect each other. Of course, once we allow 7 or more intersection 
circles to arrange themselves in 3 dimensional space, then they can mutually intersect at one point. 
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Figure 13. A Closed 5 Ring Constrained to Surround a Dimer 



Thus, so far we have established that 6 or more intersection circles can mutually intersect at no 
more than 1 point. 13 intersection circles can not mutually intersect - this is the kissing number 
limit. Thus, 2-5 intersection circles can mutually intersect at 2 points, 6-12 intersection circles 
can mutually intersect at 1 point, and > 12 intersection circles can not mutually intersect. 

In conclusion, this rule states that any adjacency matrix implying that > 6 intersection circles 
mutually intersect at 2 points is unphysical. 



This rule can be identified by the following adjacency matrix pattern: 

3 2m for which A m j = 1 for > 6 j where > 6 adjacency matrix elements amongst the j '• = 1. For 
example, j = i, l,p, q, z, y and An = Ai p = A pq = A qz = A zy = A y i = 1 (i.e. there are 6 j and a 6 
ring is formed). 



Rule 14: 

This rule states that a closed 5 ring can not be confined to surround a dimer. It is a more 
general form of the rule stating that the 5 points of a closed 5 ring can not all touch the same 
dimer (rule 9). The difference is that this rule does not require that all 5 points be touching the 
dimer, simply that it make enough contacts with the dimer such that it is confined to surround 
it; therefore forcing one distance to be < R, and thus implying particle overlap (see fig 13). 



This rule can be identified by the following adjacency matrix pattern: 
A violation occurs if 

Am = Aip = A pq = A qm = Amk = 1 (closed 5 ring) 

and 

Aij = 1 (the dimer) 

and one point in the dimer forms 3 contacts with the closed 5 ring, where 2 of the contacts are 
adjacent: Ai q = Aik = An = 1. And the other point in the dimer (here j) touches the 2 points 
of the 5 ring not touched by the other point (here i), Aj m = Aj P = 1; or touches one point not 
touched by i, and touches one of the adjacent contacts i makes, Aj m = Aji = 1 or Ajk = Aj P = 1. 



(See figure 13). If these conditions are satisfied, then the 5 ring is constrained to surround the 



dimer, and the conformation is unphysical as it implies a distance < R. 



Rule 15: 
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Figure 14. 2 Points Constrained to Opposite Sides of a Closed 4 Ring 



2 points confined to opposite sides of a closed 4 ring can not touch (see figure 14). 

This rule can be identified by the following adjacency matrix pattern: 

If A k i = An = Aij = Ajk = 1 (a closed 4 ring) 
and Akm = A m i = Aj p = A p i = 1 (2 points, m and p, are constrained to opposite sides of the 4 ring) 

T^mp > R (i-e. A mp 7^ 1 - the particles can not touch) 
If A mp = 1, this adjacency matrix is unphysical. 

Rule 16: 



This rule calculates the non-touching distances between one point added to an octahedron (See 
figure 15). This is the distance, for example, between points % and j, r^. Particles i, q, and j all 



lie on the p,k intersection circle, I p k- Thus, we can calculate r^- by determining the angle between 
Tim and Tjm, where m is the center of I p k- We can determine this angle, 6, simply by summing 
6*i and 9 2 . (We can simply sum the angles here because all of the relevent line segments, r im , 
r qmi r jm, and r^, necessarily lie in the same plane, as they all lie on the plane defined by the p, k 
intersection circle). 

#i is, as before, given by 

0i = cos" 1 Q 

(Note that this is the dihedral angle of a tetrahedron). 82 is the dihedral angle of an octahedron, 
and is given by 

1' 



6 2 = cos 1 
Thus, 6 is given by 



3 



e 1 + e 2 = cos 1 Q j +cos 1 (-- 



7T 
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Figure 15. A Particle Added to an Octahedron 
Thus, r im and r m j form a straight line, and is simply the sum of these 2 line segments: 
(23) Tij = r im + r mj = ^-R + ^-R = y/3R 

Note that, due to the symmetry of this structure, the distance between % and any point that % 
does not touch is = v3 (i.e. = the distance between i and both of the unlabeled points in the 
octahedron also). 



This rule can be identified by the following adjacency matrix pattern: 

The same pattern as in rule 6 (keeping in mind that the dummy variables used there do not 



correspond to the dummy variables shown in fig 15), and 3 1 particle, i, that touches any trimer 



within the octahedron - for example Ai q = Aik = Ai P = 1 (where q } p } k can be identified as a 
trimer because A qp = A g k = A p k = 1). In this case V iz = a/3-R, where z are the 3 points that are 



not part of the trimer to which i is attached, i.e. z = j and the other 2 unmarked points in figure 15 
Rule 17: 

The long distance, r^, between a particle added to the side of 3 connected tetrahedron (see 



figure 16) can be derived via spherical trigonometry. A spherical triangle is shown in figure 17 
Just as in standard trigonometry, there are spherical analogues to the law of cosines and the law 
of sines. In the figure, we can see that triangles OAC, OAB, and OBC form the spherical triangle 
ABC. a is the arc length opposite the angle a' (i.e. the arc length associated with the side BC), b 
is the arc opposite the angle b' (i.e. associated with the side AC, and c is the arc length opposite 
d (associated with the side AB). a' is the angle BOC (opposite vertice A), b' is the angle AOC 
(opposite vertice B), and d is the angle BOA (opposite vertice C). The arc length, S, is by 
definition given by 

S = r6 
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Figure 16. A Particle Added to the Side of 3 Connected Tetrahedron 




Figure 17. A Spherical Triangle 



where r is radius of the sphere and 9 is the angle swept out by the arc length. Therefore 

a = Ra! 
b = Kb' 
c = Rc 

where R is the radius of the sphere, and O is the center of the sphere. Thus, R is equal to the line 
segments OA, OB, and OC. (Note that, by convention, R is often taken to = 1 without loss of 
generality, which sets the arc length equal to the angle - so that a = a' and so on...). A, B, and 
C are taken both to be the vertices of the spherical triangle, as well as the dihedral angle between 
the relevant triangles, so that 

A is the dihedral angle between triangles OAC & OAB {AO AC and AOAB) 

B is the dihedral angle between AOAB & AOBC 

C is the dihedral angle between AOBC & AO AC 
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Figure 18. The Spherical Triangles Associated with the Conformation in Figure 16 



The conformation we are considering here forms 2 adjacent spherical triangles (see figure 18). 
We know the following arc lengths of the spherical triangles and their associated line segments: 
iq, qk, ik, jk. The only line segment we don't know is ij, for which we aim to solve. We also 
know the following angles (Z): Zipq, Zqpk, Zkpj. The only angle we do not know is Z.ipj, which 
we will need in order to solve for rjj. 

We can use the spherical law of cosines to calculate the associated dihedral angles. We will use 
the following form of the spherical law of cosines 

(24) sin h' sin c cos A = cos a — cos b' cos c 

Ultimately, we need to calculate the dihedral angle between Aipk and Ajpk, as this will allow 
us to back out what Z.ipj is, which we can then use to calculate the unknown distance r^. In 
spherical Aikj, the dihedral angle between between Aipk and Ajpk, Zipj, and ry are unknown. 

Thus 



This corresponds to both A and a' being unknown in the spherical rule of cosines (eqn 24). 
we can not immediately apply this equation to solve for r^. We can however use spherical Ajqk 
to calculate the dihedral angle between Aipk and Ajpk, as this dihedral angle will = the dihedral 
angle between Aipk and Akpq + the dihedral angle between Akpq and Ajpk, and spherical 
Ajqk will allow us to calculate the dihedral angle between Akpq and Ajpk. (We already know 
the dihedral angle between Aipk and Akpq; this is just twice the dihedral angle of a tetrahedron). 

Thus, let us consider spherical Ajqk; here the vertices of the spherical triangle (fig 17) corre- 
spond to the following points, A = k,B = j,C = q. Thus, A is the dihedral angle between Ajpk 
and Akpq, a' is Zjpq, b' is Zkpq, and d is Zjpk. We know all of these values except for A, they 
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arc: 



b' 



Zjpq 
Zkpq 
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3 

7T 



Zjpk = cos 



the first two angles are tt/3 as they are associated with equilateral triangles, and the last angle 
is derived from the isosceles Ajpk, where = (5/3)R has been calculated earlier in rule 4, and 



pk 



R. So, from the spherical law of cosines for spherical Ajqk, we have 



sin 



7T 



sin cos 



-7 



~2 18~~ 
5^33 
36 



cos A 

cos A 
cos A 
A 

A 



cos 

1 

2 ~ 
1 

2 ~ 

cos 



7T 

3. 

1 -7 

2 T8 
-7 
36 



cos 



7T 



COS COS 



-7 



-i 



-i 



25-36 



cos 
.514806 



36 • 5V33 
5 



This then implies that the dihedral angle between Aipk and Ajpk, <p, is given by 



2 cos 



+ cos 



where 2 cos -1 (1/3) is the dihedral angle between Akpq and Aipk (which is double the dihedral 
angle of a tetrahedron). Now, we can calcualte Zipj from the spherical law of cosines for the 
spherical Aikj, where the vertices of the spherical triangle (fig 17) correspond to the following 
points, A = k, B = j, C = i. Thus, 

b' = Zipk = — 



Zjpk 



cos 



-7 



A 



dihedral angle between Aipk and Ajpk 



2 cos 



-i 



+ cos 



-1 



7 33, 
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sin 



and a' = Zipj is unknown, and we will solve for it: 

[ — \ sin ( cos^ 1 — | cos f 2 cos -1 ( - | + cos -1 ( — = | | = cos a — cos ( — ^ cos ( cos -1 — ) 

V3 f / _i / 5 \\ ,7 

— cos 2 cos + cos —= = cos a H 

2 I 18 ) V \3J \V33j) 36 



5V33 / /1\ / 5 

cos 2 cos - + cos = cos a + 



36 V \3/ VV33// 36 

- ^ (-1) 

ps 2.94884 

and so, now that we have Zipj, we can use the regular law of cosines to calculate r^: 



Tij = K I i?2 + R2 _ 2j R2 cog ( CQS -1 



53 

54 



= i2R 2 -2R 2 



53 
54 



53 „ 



(25) 



This rule can be identified by the following adjacency matrix pattern: 

Aij = 1, and 3 4 points that touch both i and j; A P i = A P j = A q i = A q j = An = Aij = A m i = 
A m j = 1. Where 3 consecutive contacts are made amongst these 4 points; A pq = A q i = Ai m = 1 
(i.e. the same pattern as in rule 4). And 3 another particle that touches the top or bottom of the 
left or right faces (i.e. that touches either trimer Urn, ipq, jlm, or jpq. Take, for example, the 

particle z, where Ai Z = A\ z = A mz = 1, then this rule identifies V zp = |y / ^|^i?. 

Note, this packing is an enantiomer: i.e. there exist left and right handed configurations of this 
structure, which are different via chiral symmetry, but are not different in terms of their distance 
matrices. 

Rule 18: 



Following from rule 17, any adjacency matrix that implies = R is unphysical. 

This rule can be identified by the following adjacency matrix pattern: 
The same pattern as in rule 17, but where A zp = 1 also. 



natalie arkus, vinothan n. manoharan and michael p. brenner 

3. Some Comments (Non-Uniqueness of Rules) 

Rules 1-7 are sufficient to solve for packings of n < 6. For n = 7, we also need rules 
8-18. And in fact, only rules 1, 2, 4, and 6 were actually used to calculate the 6 particle 
packings (as the application of these rules alone were sufficient to either completely solve 
or eliminate as unphysical all adjacency matrices). And for 7 particles, only rules 1, 2, 3, 
4, 5, 6, 8, 9, 11, 12, 14, 16, and 17 were used (where some combination of rules 3, 5, and 9 
were used - as the way we've implemented this in the computer code, it lumps these rules 
together in the output as 'matrices were unphysical because they implied that there were 
> 2 points on an intersection circle that all touched each other'). 

Rule 17 calculates the new distance generated by adding 1 particle to a packing of 6 
particles that has already been solved. From derivations of this nature, it seems that it 
should be possible (and it would certainly be advantageous) to derive 1 general rule that 
solves for all unknown distances - at least in the most simple case of adding 1 particle to a 
packing of n particles that has already been solved. The triangular bipyramid rule derived 
in the paper is a general rule that can solve all iterative adjacency matrices (where an 
iterative matrix is one that contains only already encountered sub-matrices, i.e. contains 
no new structure). Thus, all iterative rules derived up until this point can be replaced by 
this rule - this is one example of the non-uniqueness of rules that can be used to derive 
packings; however, the non-iterative rules (those applied to physical as well as unphysical 
new seeds) must still be applied. 

For n > 8 the triangular bipyramid rule derived in the paper was used to solve for all it- 
erative adjacency matrices and individual geometrical rules were derived for non-iterative 
adjacency matrices. The following table shows that iterative matrices account for most of 
the adjacency matrices, and that the number of non-iterative adjacency matrices does not 
become 'too large' until n = 10. 



Number of 


Number of 


Iterative 


Non-Iterative 


Particles 


Adjacency Matrices 


Adjacency Matrices 


Adjacency Matrices 


6 


4 


3 


1 


7 


29 


26 


3 


8 


438 


437 


1 


9 


13,828 


13,823 


5 


10 


750,352 


750,258 


94 



At n = 10, there are 94 non-iterative adjacency matrices. We have solved for the packings 
that correspond to these matrices numerically, which is why the new seed list at n = 10 is 
only a preliminary list. Future work will derive a complete list of n = 10 new seeds, as well 
as go up to higher n. One way of accomplishing this is to apply the triangular bipyramid 
rule to new seeds. 

New seed rules for n = 8 and 9 can be derived similarly. Please refer to the computer code 
for the relevant rules. The code can be downloaded from 

http : //people . seas . harvard . edu/~narkus/pubs . html 
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and is entitled 'EvalAdjMats.c.' In the code, solutions for 8 and 9 particle new seeds were 
input from Newton iterations to a precision of 16 decimal places. 

4. Derivation of the Triangular Bipyramid Rule 

We derive this rule as follows: In 3 dimensions, let the distance, r^, between 2 points, i and j, 
be unknown. The points i and j share a common 3 particle base, p, k, q. The 5 points i, j, k,p, q 



together form a (potentially irregular) ditetrahedron or triangular bipyramid (see figure 19). This 



triangular bipyramid can be decomposed into 3 related tetrahedra (see figure 20). Let us first 



derive a formula that relates these dihedral angles (Ai, Bi, C{) to the angles (a;, bi, q); I — 1, 2, 3. 



Consider the general (potentially irregular) tetrahedron depicted in figure 21 The dihedral 
angle, A, between planes AOB and AOC can be calculated using the dot product of the normals 
to the planes. The normals to the planes are given by the cross products of the vectors to the 
vertices. Thus, this dot product is given by 

(OA x OB) ■ (OA x OC) = (\OA\\OB\smc)(\OA\\OC\smb)cosA 

(26) = | OA\ 2 1 OB 1 1 OC\ sine sin b cos A. 
Using a well known vector identity, we also know that 

(OA x OB) ■ (OA x OC) = OA ■ [OB x (OA x OC)] 

= OA ■ [OA(OB ■ OC) - OC(OA ■ OB)] 

= (OB ■ OC) — (OA ■ OC)(OA ■ OB) 

= \OB\\OC\ cosa - \OA\\OC\ cosb\OA\\OB\ cose 

(27) = \OB\\OC\ cos a -\OA\ 2 \OC\\OB\cosb cose 



Setting eqn 26 = eqn 27, we thus have 



|OA| 2 |OS||OC|sincsin6cosA = \OB\\OC\ cosa - \OA\ 2 \OC\\OB\ cosfe cose 
We can divide by \OB\ and \OC\, leaving us with 

|Ov4| 2 sin c sin 6 cos A = cos a — \OA\ 2 cos b cose 
and without loss of generality, we can always rescale things so that \OA\ = 1, leaving us with 

(28) sin c sin b cos A = cos a — cos b cos c 

Note that this formula is simply the spherical rule of cosines; however, in this derivation we do 
not assume that all radial distances are equal, as is customary in the derivation of the spherical 
rule of cosines. This can analogously be derived for dihedral angles B and C, in which case we 
obtain 

(29) sin c sin a cos B = cos b — cos a cos c 

(30) sin a sin b cos C = cos c — cos a cos b 



To derive the general formula for r^, we begin with the 3rd tetrahedron of Fig. 20 (that with 



angles A 3 , B 3 , a 3 , 6 3 , etc.) and write down the expression for A 3 using eqn 28 

_ x f cos a 3 — cos b 3 cos c 3 



A 3 = cos 



sin c 3 sin b 3 
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k 



q 



Figure 19. The Triangular Bipyramid. 

The triangular bipyramid (or ditetrahedron) constructed in the triangular bipyramid 
rule. The center triangle (kpq), shown in red, corresponds to the common 3 particle 
base. Particles % and j are related to one another through the common base. The 
distance between i and j, r^, shown as the dash-dot blue line, is unknown. a% 
corresponds to ^jpi. A\ is the dihedral angle between Aipk and Ajpk, A 2 is the 
dihedral angle between Aipk and Akpq, and A3 is the dihedral angle between Akpq 
and Ajpk. Points i and j can either both lie on the same side of the base kpq or 
each lie on opposite sides of the base (indicated by the dashed lines, that can either 
go into or come out of the plane). If i,j lie on the same side, then A 1 is equal to 
the difference of A 2 and A 3 , and if i,j lie on opposite sides of the base then Ai is 
equal to the sum of A2 and A3. When all distances other than are known, then 
an analytical formula can be derived to solve r^. 

Analogously, for the 2nd tetrahedron, we have 

1 / cos a 2 — cos b 2 cos c 2 \ 



A\ will be either the sum or the difference of the dihedral angles A 2 and A 3 (see fig [19] ), 
depending on whether the points i,j lie on the same or on opposite sides of the base p, k, q (on 
the same side, A\ corresponds to the difference, and on opposite sides to the sum). 

a,\ is then given by 




(31) ai = cos (sinci sin 61 cos A\ + cos 61 cosci) 

and finally, from the law of cosines, we can calculate r^: 
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Figure 20. The Tetrahedra for the Triangular Bipyramid Rule. The 3 

tetrahedra constructed for the triangular bipyramid rule. When applied to iterative 
packings, all of the sides of the tetrahedra are known except for side ij (correspond- 
ing to the distance r^). The known distances can be used to calculate A 2 and A 3 , 
which in turn allows us to calculate A±, and thus to solve for r\y 
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Figure 21. A General Tetrahedron. 

A, .B,C and O label the points of the tetrahedron. A is also used to correspond 
to the dihedral angle between AAOB and AAOC, B corresponds to the dihedral 
angle between AAOB and ABOC, and C corresponds to the dihedral angle between 
ABOC and AAOC. a corresponds to ZCOA, b to ZCOB, and c to AAOB. 

Associated with each we have 2 possible A ± , and thus 2 possible solutions (similar, in 
principle, to one having 2 possible solutions to a quadratic equation). 

4.1. Derivation of the Formula To be Used for New Seeds. From setting Ai either equal 
to the sum or the difference of A 2 and A 3 , we have 



cos 



cos ai — cos bi cos ci 



= cos 



-l 



cos a 2 — cos 6 2 cos c 2 



± cos" 



cos a 3 — cos 6 3 cos c 3 



sin c 2 sin b 2 J \ sin c 2 sin b 2 J \ sin c 3 sin 6 3 

writing the angles a iy foj, q in terms of their relevant distances, using the law of cosines, we then 
have 



cos 



d\+d\-d\ 
2d^d4 



2didfj 2d\d^ 



_ (dl+dl-d^ 
Mfdj 



1 Ad\d\ 



COS 
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and solving for d 2 then yields 



— 2rf 6 C?4A / 1 — 
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APPENDIX B: PACKINGS OF 8 < n < 10 SPHERES 



Here we include the complete set of packings for 8 < n < 10 particles, derived via the method 
specified in the paper. (Packings of n < 7 particles can be found in the paper). Throughout, R 
corresponds to twice the radius of the particles, and without loss of generality can be normalized 
to 1, if desired. In all instances, a '*' appears in front of the 2nd moment that corresponds to 
the minimum of the 2nd moment of n particles. If a packing is denoted as chiral, then it has 2 
distinct states. 

Note that the code used to generate these results as well as the raw output from the code (in 
which case the results are output to 16 decimal places instead of the 4 d.p. shown here) can be 
downloaded from the following website: http://people.seas.harvard.edu/~narkus/. 

Notation: 

A: Adjacency Matrix. 

T>: Distance Matrix. 

C: Coordinates of a packing. 

0: Point group of a packing. 

a: Symmetry number of a packing. 

1. Packings of n = 8 Particles 

Here, we include the 13 packings of n = 8 particles. The 'Special Properties' column denotes 
whether the packing is chiral or corresponds to a new seed. If the special properties column is 
blank, then the packing is neither. 1 of these packings (packing 13) corresponds to a new seed, 
and 3 of the packings are chiral. The minimum of the 2nd moment corresponds to packing 13, 
graph 408. 
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Packing 1 (Graph 27): 
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Packing 2 (Graph 77): 
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Packing 3 (Graph 82): 



V : 



(° 








1 





1 


1 


1\ 














1 


1 





1 














1 





1 


1 


1 

















1 


1 





1 


1 








1 


1 


1 


1 


1 








1 





1 


1 


1 





1 


1 


1 


1 





1 


\1 


1 


1 


1 


1 


1 


1 





( 


1 


6667 


1.6667 


1 


1.6330 


1 


1 


A 


1.6667 







1.6667 


1.9907 


1 


1 


1.6330 


l 


1.6667 


1 


6667 





1.0887 


1 


1.6330 


1 


i 


1 


1 


9907 


1.0887 





1.6667 


1.6330 


1 


i 


1.6330 




1 


1 


1.6667 





1 


1 


l 


1 




1 


1.6330 


1.6330 


1 





1 


l 


1 


1 


6330 


1 


1 


1 


1 





l 


V i 




1 


1 


1 


1 


1 


1 


0/ 



R 








1.6667 


1.4434 
0.8333 


0.9302 
-0.0556 
-0.3629 
0.7698 
1.3333 
0.5443 
-0.0962 
0.8333 
0.5443 
0.7698 
0.3333 
0.5443 
0.4811 
0.8333 
V-0.2722/ 



R 




2nd Moment 





a 


Special Properties 


5.61574R 2 


Ci 


1 


chiral 



Packing 4 (Graph 107): 
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Packing 5 (Graph 204): 
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Packing 6 (Graph 224): 
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Packing 7 (Graph 271): 
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Packing 8 (Graph 294): 
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Packing 9 (Graph 358): 
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Packing 10 (Graph 366): 
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Packing 11 (Graph 374): 
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Packing 12 (Graph 397): 
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Packing 13 (Graph 408): 
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2. Packings of n = 9 Particles 

Here, we include the 52 packings of n = 9 particles. The 'Special Properties' column denotes 
whether the packing is chiral, non-rigid, or corresponds to a new seed. If the special properties 
column is blank, then the packing is none of these. 28 of these packings are chiral, 4 of these pack- 
ings correspond to new seeds (graphs 8901, 8926, 12811, and 13380), and 1 of them is non-rigid 
(graph 8901). The minimum of the 2nd moment corresponds to packing 47, graph 13380. A star, 
'*', appears in front of packings that are reported here but were not reported in our last paper 
(PRL, 103, 118303). 



Packing 1 (Graph 406): 
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Packing 2 (Graph 756): 
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Packing 3 (Graph 971): 
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Packing 4 (Graph 991): 
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Packing 5 (Graph 1039): 
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Packing 6 (Graph 1050): 
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Packing 7 (Graph 1427): 
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Packing 8 (Graph 1946): 
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Packing 10 (Graph 3313): 
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Packing 12 (Graph 3850): 
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Packing 14 (Graph 4549): 
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Packing 15 (Graph 4564): 
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Packing 16 (Graph 4574): 



A : 



V : 











1 





1 





1 
















1 





1 


1 


1 

















1 


1 


1 


1 


1 














1 





1 








1 














1 





1 


1 





1 


1 











1 


1 





1 


1 





1 








1 


1 


1 


1 


1 


1 





1 


1 





1 


\1 


1 


1 





1 


1 


1 


1 


o/ 



/ 

1.0887 
1.6330 
1 

1.7249 
1 

1.6667 
1 

V i 



1.0887 


1.6330 
1.7249 
1 

1.6667 
1 
1 
1 



1.6330 
1.6330 


1.6667 
1.6667 

1 

1 

1 

1 



1 

1.7249 
1.6667 


2.4637 
1 

1.9907 
1 

1.6330 



1.7249 

1 

1.6667 
2.4637 


1.9907 
1 

1.6330 
1 



1 

1.6667 
1 
1 

1.9907 


1.6330 
1 
1 



1.6667 
1 
1 

1.9907 
1 

1.6330 

1 
1 



1 
1 
1 
1 

1.6330 
1 
1 

1 



1 - 

1 

1 

1.6330 
1 
1 
1 
1 

, 



R 




2nd Moment 





a 


Special Properties 


7.39735.R 2 


c 2 


2 


chiral 



DERIVING FINITE SPHERE PACKINGS 



103 



Packing 17 (Graph 4811): 
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Packing 18 (Graph 4816): 
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Packing 19 (Graph 5093): 
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Packing 20 (Graph 5098): 
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Packing 21 (Graph 5101): 
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Packing 22 (Graph 5112): 
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Packing 23 (Graph 5116): 
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Packing 24 (Graph 5399): 



A : 



V 











1 


1 





1 



















1 


1 


1 


1 


1 

















1 





1 


1 


1 

















1 





1 


1 


1 














1 


1 


1 





1 


1 














1 


1 


1 


1 





1 


1 











1 





1 


1 





1 


1 








1 


\1 


1 


1 


1 


1 


1 


1 


1 


0/ 



/ 

1.6330 
1.7778 
1 

1 

1.9907 
1 

1.6667 

V i 



1.6330 


1.6330 
1.6667 

1 

1 

1 

1 

1 



1.7778 
1.6330 


1.8178 
1.6667 
1 

1.9907 
1 
1 



1 

1.6667 
1.8178 


1.6330 
1.7778 
1 

1.9907 
1 



1 
1 

1.6667 
1.6330 


1.6330 
1 
1 
1 



1.9907 
1 
1 

1.7778 
1.6330 


1.6667 
1 
1 



1 
1 

1.9907 
1 
1 

1.6667 


1.6330 
1 



1.6667 
1 
1 

1.9907 
1 
1 

1.6330 

1 



1 
1 
1 
1 
1 
1 
1 

0/ 



R 



C : 








1.6330 


-1.4913 
0.9677 


-0.1104 
0.2722 
-0.9559 


0.8165 
0.5735 

-0.9345 
1.7237 

-0.3441 
0.4636 
0.8165 

-0.3441 

-0.7726 
1.3608 
0.5735 

-0.5298 
0.8165 

-0.2294/ 



R 




2nd Moment 





a 


Special Properties 


7.43987.R 2 


c 2 


2 


chiral 



DERIVING FINITE SPHERE PACKINGS 



111 



*Packing 25 (Graph 5416): 
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Packing 26 (Graph 5544): 
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Packing 27 (Graph 5835): 
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Packing 28 (Graph 6031): 
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Packing 29 (Graph 8901): 
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Packing 30 (Graph 8926): 
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Packing 31 (Graph 8943): 
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Packing 32 (Graph 8955): 
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Packing 33 (Graph 8959): 
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Packing 34 (Graph 9352): 
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Packing 35 (Graph 9503): 
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Packing 36 (Graph 9520): 
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Packing 37 (Graph 10279): 
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*Packing 38 (Graph 10664): 
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Packing 39 (Graph 10989): 
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Packing 40 (Graph 10997): 
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Packing 41 (Graph 11229): 
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Packing 42 (Graph 11245): 
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Packing 43 (Graph 11633): 
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Packing 44 (Graph 11998): 
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Packing 45 (Graph 12163): 
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Packing 46 (Graph 12180): 
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Packing 47 (Graph 12239): 
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Packing 48 (Graph 12811): 
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Packing 49 (Graph 13380): 
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Packing 50 (Graph 13521): 
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Packing 51 (Graph 13576): 
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Packing 52 (Graph 13700): 
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3. Packings of n 



10 Particles 




Here we include the 262 10 particle packing^ \ The 'special properties' column denotes whether 
a packing is chiral, corresponds to a new seed, to a non-rigid structure, or to a packing with > 3n— 6 
contacts. If this column is blank, then that packing has none of these properties. There are 201 
chiral packings, 8 new seeds (graphs 308923, 618718, 622587, 629072, 665570, 679397, 714961, and 
749808), 4 non-rigid packings (Graph 665570 corresponds to the packing that is both a new seed 
and a non-rigid structure, and graphs 286981, 457749, and 593702 (packings 84, 145, and 197) 
correspond to the 10 particle non-rigid packings that are derived from the 9 particle non-rigid new 
seed; 9 particle graph 89010 and 3 packings with 3n — 5 = 25 contacts (packings 259-262). The 
minimum of the 2nd moment corresponds to packing 257, graph 749808. 

Note that for the 25 contact packings, we have included all graphs that lead to these packings; 
thus multiple graphs are listed under the same packing header. Because of this, we have listed all 
25 contact packings at the end; these are the only packings that are listed out of chronological 
graph order. For these graphs, we have included their associated distance matrix and coordinates, 
as these can be different up to particle labeling degeneracy (i.e. they can be isomorphic but not 
identical), but we have only included the table with the 2nd moment, 0, cr, and the special 
properties once, as this will always be identical. 

A star, '*', appears in front of packings that are reported here but were not reported in our last 
paper (PRL, 103, 118303). A double star, '**', appears in front of packings that are reported here 
but were not reported previously by us or by Hoy and O'Hern (PRL, 105, 068001). We have also 
included markings next to the relevant graphs that correspond to 25 contact packings which were 
not found in these papers; however, note that this did not affect the number of packings reported 
in either of these papers as multiple graphs lead to the same 25 contact packings, and in both 
papers all 3 25 contact packings were found. In front of these graphs, a '*' appears for ones not 



found in (PRL, 103, 118303), and a '+' appears for ones not found in (PRL, 105, 0680010 



Note that for these 10 particle packings, the packing figures were generated by code that John Lee wrote that 
interfaces the packing output with the mayavi2 graphing package. And the point groups and symmetry numbers 
were also obtained via code that John Lee wrote that interfaces the packing output with the rotational constant 
calculator website (http://www.colby.edu/chemistry/PChem/scripts/ABC.html) that calculates these numbers. 

26 Note that while the 3 non-rigid packings based on the 9 particle non-rigid new seed were reported in (PRL, 
103, 118303) and the PhD Thesis (http://people.seas.harvard.edu/~narkus/assets/Thesis.pdf.zip), they were not 
found by the old code. We had reported them because we had already solved for their coordinates. We made a note 
of this where the packings were listed, and mentioned that we suspected these 3 packings had been missed due to 
round-off error, which we knew to be an issue. Thus, while a star does not appear next to these packings because 
they were previously reported, if one were to compare raw output of the old and new code, these 3 matrices would 
also lie in the category that the old code misses and the new modified code picks up. 

27 Rob Hoy sent us the raw data he used for the paper (PRL, 105, 068001). It was comparing this data to our 
data that allowed us to determine which graphs (ones that uniquely or redundantly correspond to packings) appear 
here but were not reported in (PRL, 105, 068001). 
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Packing 1 (Graph 4471): 
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Packing 3 (Graph 17438): 
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Packing 4 (Graph 18823): 
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Packing 5 (Graph 23979): 
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Packing 6 (Graph 26619): 
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Packing 7 (Graph 26632): 
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Packing 8 (Graph 29760): 
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Packing 9 (Graph 31998): 
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Packing 10 (Graph 32015): 
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Packing 11 (Graph 34257): 
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Packing 12 (Graph 40878): 
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*Packing 13 (Graph 53546): 
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Packing 14 (Graph 57812): 
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Packing 15 (Graph 62297): 
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Packing 16 (Graph 65350): 
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Packing 17 (Graph 65352): 
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Packing 18 (Graph 65490): 
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Packing 19 (Graph 69229): 
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Packing 20 (Graph 69305): 
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Packing 21 (Graph 75933): 
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Packing 22 (Graph 76012): 
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Packing 23 (Graph 78584): 
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Packing 24 (Graph 78598): 
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Packing 25 (Graph 78941): 



A : 



V : 



/0 








1 





1 





1 




























1 





1 


1 


1 


























1 





1 





1 





























1 


1 


1 


1 














1 














1 





1 

















1 


1 














1 





1 














1 





1 


1 











1 


1 














1 


1 


1 





1 








1 


1 














1 1 





1 


1 





1 


1 





1 














Vi i 


1 


1 





1 


1 


1 


1 


o) 














/ 




1.0887 


1.9959 


1.6330 




1 


1.7249 


1 


1.6667 


1 




1.0887 









1.6330 


1.6330 


1.7249 


1 


1.6667 


1 


1 





1.9959 




1.6330 







1.0887 


2.4815 


1 


1.7249 


1 


1.6667 





1.6330 




1.6330 


1.0887 







1.6667 


1.6667 


1 


1 


1 





1 




1.7249 


2.4815 


1.6667 







2.4637 


1 


1.9907 


1 





1.7249 




1 




1 




1.6667 


2.4637 





1.9907 


1 


1.6330 





1 




1.6667 


1.7249 




1 




1 


1.9907 





1.6330 


1 





1.6667 




1 




1 






1 


1.9907 


1 


1.6330 





1 





1 




1 




1.6667 




1 




1 


1.6330 


1 


1 








V 




































0/ 



R 



* 




-1.0887 


1.6319 
-1.1491 


1.2912 
-0.5443 
-0.8386 
-0.1042 

0.3629 
-0.9260 

0.8036 
-1.4515 

0.4717 

0.8070 

0.2722 
-0.5241 

0.8070 
-1.3608 
-0.5241 

0.2925 
-0.5443 
-0.7862 

0.8372 
-0.5443 
, 



R 




2nd Moment 





a 


Special Properties 


8.93224R 2 




1 


chiral 



DERIVING FINITE SPHERE PACKINGS 



165 



Packing 26 (Graph 78955): 
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Packing 27 (Graph 80933): 
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Packing 29 (Graph 81182): 
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Packing 30 (Graph 82672): 
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Packing 31 (Graph 93924): 
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Packing 32 (Graph 93930): 
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Packing 33 (Graph 93942): 
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Packing 34 (Graph 94152): 
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*Packing 35 (Graph 94932): 
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Packing 37 (Graph 113048): 
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Packing 39 (Graph 119213): 
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*Packing 40 (Graph 138747): 
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Packing 41 (Graph 165946): 
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Packing 43 (Graph 184811): 
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*Packing 44 (Graph 187508): 
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Packing 45 (Graph 200848): 
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Packing 47 (Graph 200875): 
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Packing 48 (Graph 200880): 
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Packing 53 (Graph 201576): 
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Packing 55 (Graph 203149): 
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Packing 57 (Graph 203204): 
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Packing 58 (Graph 203235): 
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Packing 59 (Graph 203248): 
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Packing 61 (Graph 203255): 
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Packing 62 (Graph 203262): 
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Packing 63 (Graph 204149): 
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Packing 64 (Graph 204201): 
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Packing 65 (Graph 209891): 
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*Packing 66 (Graph 238324): 
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*Packing 67 (Graph 240312): 
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Packing 68 (Graph 266450): 
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Packing 69 (Graph 272310): 
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Packing 70 (Graph 282787): 
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Packing 71 (Graph 282820): 
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Packing 73 (Graph 283029): 
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Packing 74 (Graph 283324): 
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Packing 75 (Graph 283448): 
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Packing 76 (Graph 283463): 
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Packing 77 (Graph 285359): 
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Packing 78 (Graph 285374): 
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Packing 79 (Graph 286096): 
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Packing 80 (Graph 286101): 
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Packing 81 (Graph 286126): 
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Packing 83 (Graph 286160): 



A : 



V : 



/0 





1 





1 








!\ 




















1 








1 


1 


1 























1 


1 





1 


























1 





1 





1 














1 1 











1 





1 


1 

















1 


1 








1 


1 


1 














1 


1 


1 











1 


1 














1 





1 


1 








1 


1 














1 


1 


1 


1 


1 


1 





1 














Vi i 





1 1 


1 


1 


1 


1 


o) 














/ 




1.6667 


1.9868 


1.8144 




1 


1.7249 


1 


1.9907 


1.6330 




1.6667 







1.9868 


1.6667 




1 


1.6330 


1.6330 


1 


1 





1.9868 




1.9868 







1.9868 


1.6859 


1 


1 


1.6859 


1 





1.8144 




1.6667 


1.9868 







1.9907 


1 


1.7249 


1 


1.6330 





1 




1 


1.6859 


1.9907 







1.6667 


1 


1.6330 


1 





1.7249 




1.6330 


1 






1 


1.6667 





1.0887 


1 


1 





1 




1.6330 


1 




1.7249 




1 


1.0887 





1.6667 


1 





1.9907 




1 


1.6859 




1 


1.6330 


1 


1.6667 





1 





1.6330 




1 


1 




1.6330 




1 


1 


1 


1 








V 


































0/ 



R 








1.6667 


1.8036 
0.8333 


0.4563 
0.9877 
1.4521 
0.1986 
0.8333 

-0.5159 
1.2141 
0.9259 
0.8025 
0.9403 
0.3333 

-0.0688 
0.6093 
1.7222 
0.7910 
0.9403 
1.3333 

-0.0688 
0.2716 
0.8333 
V 0.4815 / 



R 




2nd Moment 





a 


Special Properties 


8.69036/? 2 


C s 


1 





DERIVING FINITE SPHERE PACKINGS 



223 



Packing 84 (Graph 286981): 
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Packing 85 (Graph 287098): 
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Packing 86 (Graph 287190): 
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Packing 87 (Graph 287227): 
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Packing 88 (Graph 287262): 
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Packing 89 (Graph 287270): 
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Packing 90 (Graph 287539): 
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*Packing 91 (Graph 287604): 
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*Packing 93 (Graph 287714): 
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Packing 94 (Graph 287739): 
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Packing 95 (Graph 290619): 
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Packing 96 (Graph 291717): 
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Packing 97 (Graph 291725): 



A : 



V : 



/0 





1 





1 








!\ 




















1 








1 


1 


























1 


1 





1 


1 























1 





1 





1 














1 1 











1 


1 


1 


1 

















1 


1 








1 


1 


1 














1 


1 


1 











1 


1 














1 





1 1 


1 








1 


1 














1 


1 


1 


1 


1 


1 




















Vi o 


1 


1 1 


1 


1 


1 





o) 














/ 




1.9938 


1.6779 


1.8155 




1 


1.9843 


1 


1.6779 


1.6532 




1.9938 







1.9843 


1.9938 




1 


1.6779 


1.6779 


1 


1 





1.6779 




1.9843 







1.6779 


1.6180 


1 


1 


1.6180 


1 





1.8155 




1.9938 


1.6779 







1.6779 


1 


1.9843 


1 


1.6532 





1 




1 


1.6180 


1.6779 







1.6180 


1 


1 


1 





1.9843 




1.6779 


1 






1 


1.6180 





1.6180 


1 


1 





1 




1.6779 


1 




1.9843 




1 


1.6180 





1.6180 


1 





1.6779 




1 


1.6180 




1 




1 


1 


1.6180 





1 





1.6532 




1 


1 




1.6532 




1 


1 


1 


1 








V 


































0/ 



R 



' 




1.9938 


-1.5177 
0.7155 


-0.6962 
0.8265 
1.4588 


0.9969 
-0.0214 
-1.2926 

1.2783 

0.7954 
-0.6722 

0.5416 
-0.5048 
-0.3079 

1.4522 

0.7822 
-0.8236 

1.4315 
-0.0744 
-0.6624 

0.5623 
V 0.4950 / 



R 




2nd Moment 





a 


Special Properties 


8.63469i? 2 


C s 


1 





DERIVING FINITE SPHERE PACKINGS 237 



Packing 98 (Graph 291912): 
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Packing 99 (Graph 291998): 
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Packing 100 (Graph 292451): 



A : 



V 



/O 





1 





1 








!\ 




















1 








1 





1 























1 


1 





1 


1 























1 





1 


1 


1 














1 1 











1 


1 





1 

















1 


1 











1 


1 














1 


1 


1 








1 


1 


1 














1 





1 1 





1 





1 


1 

















1 


1 


1 


1 


1 




















Vi i 


1 


1 1 


1 


1 


1 





o) 














/ 




1.6667 


1.1220 


1. 


9889 




1 


1.7566 


1 


1.6330 


1.7007 




1.6667 







1.9889 


1.1220 




1 


1.7566 


1.6330 


1 


1.7007 





1.1220 




1.9889 







1.6180 


1.6779 


1 


1 


1.6180 


1 





1.9889 




1.1220 


1.6180 







1.6779 


1 


1.6180 


1 


1 





1 




1 


1.6779 


1.6779 







1.9843 


1 


1 


1.6532 





1.7566 




1.7566 


1 






1 


1.9843 





1.6180 


1.6180 


1 





1 




1.6330 


1 




1.6180 




1 


1.6180 





1 


1 





1.6330 




1 


1.6180 




1 




1 


1.6180 


1 





1 





1.7007 




1.7007 


1 






1 


1.6532 


1 


1 


1 








V 


































0/ 



R 



' 




-1.6667 


1.1218 
-0.0243 


1.1218 
-1.6424 


-0.2661 
-0.8333 
-0.4845 

1.4727 
-0.8333 

0.4715 

0.5539 
-0.3333 
-0.7629 

0.5539 
-1.3333 
-0.7629 

1.3866 
-0.8333 
-0.5248 

0.5431 
-0.8333 
V 0.1030 / 



R 




2nd Moment 





a 


Special Properties 


8.09937i? 2 


C s 


1 





240 



NATALIE ARKUS, VINOTHAN N. MANOHARAN AND MICHAEL P. BRENNER 



*Packing 101 (Graph 294195): 
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*Packing 102 (Graph 294321): 
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Packing 103 (Graph 294347): 
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Packing 104 (Graph 297050): 
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*Packing 105 (Graph 298453): 
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Packing 106 (Graph 308923): 
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Packing 107 (Graph 310808): 
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Packing 108 (Graph 317168): 
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*Packing 109 (Graph 318627): 
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*Packing 110 (Graph 319641): 
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Packing 111 (Graph 364071): 
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Packing 112 (Graph 379658): 
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Packing 113 (Graph 394490): 
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Packing 114 (Graph 394589): 
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Packing 115 (Graph 396368): 
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Packing 116 (Graph 396440): 
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Packing 117 (Graph 396446): 
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Packing 118 (Graph 396466): 
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Packing 119 (Graph 396840): 
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Packing 120 (Graph 396933): 
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Packing 121 (Graph 398930): 
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Packing 122 (Graph 401472): 
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Packing 123 (Graph 420068): 
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Packing 124 (Graph 420079): 
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Packing 125 (Graph 420083): 
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Packing 126 (Graph 420106): 
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Packing 127 (Graph 420160): 
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Packing 128 (Graph 420180): 
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Packing 129 (Graph 420187): 
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Packing 130 (Graph 420198): 
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Packing 131 (Graph 422919): 
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Packing 132 (Graph 423191): 
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Packing 133 (Graph 423197): 
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Packing 134 (Graph 424068): 
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Packing 135 (Graph 424259): 
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*Packing 136 (Graph 427393): 
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*Packing 137 (Graph 445153): 
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Packing 138 (Graph 449382): 
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Packing 139 (Graph 449396): 
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Packing 140 (Graph 449426): 
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Packing 141 (Graph 449737): 
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Packing 142 (Graph 449787): 
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*Packing 143 (Graph 454720): 
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Packing 144 (Graph 456547): 
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Packing 145 (Graph 457749): 
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Packing 146 (Graph 457823): 
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Packing 147 (Graph 457894): 
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Packing 148 (Graph 458045): 
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Packing 149 (Graph 458120): 
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Packing 150 (Graph 458161): 
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Packing 151 (Graph 458198): 
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Packing 152 (Graph 459996): 



A : 



/O 





1 





1 





1 





!\ 























1 





1 


1 


1 


1 


























1 


1 





1 


1 














1 

















1 





1 














1 














1 





1 

















1 


1 














1 


1 


1 














1 


1 





1 











1 


1 














1 1 





1 





1 








1 


1 














1 


1 





1 


1 


1 


1 




















Vi i 


1 


1 





1 


1 


1 





o) 














/ 




1.6779 


1.6779 




1 


2.4520 


1 


1.9843 


1 


1.6532 




1.6779 









1.6180 


1.1220 




1 


1.6180 


1 


1 


1 





1.6779 




1.6180 







1.9889 


1.6779 


1 


1 


1.6180 


1 





1 




1.1220 


1.9889 







2.1063 


1.6330 


1.7566 


1 


1.7007 





2.4520 




1 




1.6779 


2.1063 







1.9843 


1 


1.6779 


1 





1 




1.6180 


1 




1.6330 


1.9843 





1.6180 


1 


1 





1.9843 




1 




1 




1.7566 




1 


1.6180 





1.6180 


1 





1 




1 




1.6180 




1 


1.6779 


1 


1.6180 





1 





1.6532 




1 




1 




1.7007 




1 


1 


1 


1 








V 




































0/ 



R 



' 




-1.6779 


1.4175 
-0.8978 


-0.5319 
-0.7618 

0.3697 

0.6434 
-2.3326 
-0.3969 

0.7671 
-0.3568 
-0.5331 

0.9435 
-1.7143 

0.3295 
-0.1090 
-0.8390 
-0.5331 

0.7459 
-1.3554 
-0.5827 

0.4617 
-0.8390 
V 0.2880 / 



R 




2nd Moment 





a 


Special Properties 


8.71068i? 2 




1 


chiral 



292 



NATALIE ARKUS, VINOTHAN N. MANOHARAN AND MICHAEL P. BRENNER 



*Packing 153 (Graph 460165): 
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Packing 154 (Graph 460203): 
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Packing 155 (Graph 460314): 
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*Packing 156 (Graph 460524): 
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*Packing 157 (Graph 460588): 
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Packing 158 (Graph 462814): 
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Packing 159 (Graph 463479): 
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*Packing 160 (Graph 463835): 
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Packing 161 (Graph 464807): 
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Packing 162 (Graph 466085): 
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Packing 163 (Graph 466132): 
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Packing 164 (Graph 466155): 
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Packing 165 (Graph 466323): 
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Packing 166 (Graph 466327): 



A : 



/O 





1 


1 








1 


!\ 




















1 





1 


1 





1 























1 


1 





1 


1 














1 














1 


1 

















1 











1 


1 




















1 


1 














1 


1 














1 


1 


1 








1 


1 


1 














1 





1 1 





1 





1 


1 














1 


1 


1 


1 


1 


1 





1 














Vi i 


1 





1 


1 


1 


1 


o) 














/ 


1 


7249 


1.6330 




1 


2.0718 


1 


1.6667 


1.0887 


1 




1.7249 







1.6667 


1.9868 




1 


1.9907 


1 


1 


1.6330 





1.6330 


1 


6667 







1.9868 


1.9907 


1 


1 


1.6330 


1 





1 


1 


9868 


1.9868 







1.8113 


1.6859 


1.6859 


1 


1 





2.0718 




1 


1.9907 


1.8113 







2.4369 


1 


1 


1.6667 





1 


1 


9907 


1 




1.6859 


2.4369 





1.6330 


1.6667 


1 





1.6667 




1 


1 




1.6859 




1 


1.6330 





1 


1 





1.0887 




1 


1.6330 




1 




1 


1.6667 


1 





1 





1 


1 


6330 


1 






1 


1.6667 


1 


1 


1 








V 


































0/ 



R 




2nd Moment 





a 


Special Properties 


9.01999i? 2 




1 


chiral 



306 



NATALIE ARKUS, VINOTHAN N. MANOHARAN AND MICHAEL P. BRENNER 



Packing 167 (Graph 466457): 
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*Packing 168 (Graph 466903): 



A : 



/0 





1 


1 





1 


1 


!\ 




















1 





1 





1 


1 























1 


1 


1 





1 














1 

















1 


1 














1 











1 





1 

















1 


1 











1 





1 














1 


1 


1 








1 





1 














1 


1 





1 


1 





1 


1 














1 1 





1 1 








1 





1 














Vi i 


1 


1 


1 


1 


1 


1 


0/ 














/ 




1.6859 


1.6330 




1 


1.9770 


1 


1.6667 


1 


1 




1.6859 







1.6859 


1.1471 




1 


1.9868 


1 


1.6059 


1 





1.6330 




1.6859 







1.9907 


1.9770 


1 


1 


1 


1.6667 





1 




1.1471 


1.9907 







1.7791 


1.6667 


1.7249 


1.6330 


1 





1.9770 




1 


1.9770 


1.7791 







2.3778 


1 


1.5789 


1 





1 




1.9868 


1 




1.6667 


2.3778 





1.6330 


1 


1.6330 





1.6667 




1 


1 




1.7249 




1 


1.6330 





1 


1.0887 





1 




1.6059 


1 




1.6330 


1.5789 


1 


1 





1 





1 




1 


1.6667 




1 




1 


1.6330 


1.0887 


1 








V 


































0/ 



R 



\ 






1.6859 


1.4287 
0.7909 


-0.5185 
0.7493 

-0.4120 

-0.0109 
1.7056 
0.9997 
0.9505 

-0.0312 

-0.3090 
0.7969 
1.3702 
0.5151 
0.7259 
0.3746 
0.5769 

-0.1555 
0.8429 
0.5151 
0.4666 
0.8429 

-0.2678/ 



R 




2nd Moment 





a 


Special Properties 


8.6358CLR 2 




1 


chiral 



308 



NATALIE ARKUS, VINOTHAN N. MANOHARAN AND MICHAEL P. BRENNER 



Packing 169 (Graph 466969): 
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Packing 170 (Graph 467170): 
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Packing 171 (Graph 467260): 
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Packing 172 (Graph 472149): 



A : 



f° 








1 





1 





1 


1 
















1 





1 


1 





1 

















1 


1 


1 


1 





1 























1 


1 





1 














1 








1 


1 





1 














1 


1 








1 


1 





1 








1 


1 


1 


1 


1 


1 








1 


1 








1 


1 





1 


1 





1 


1 








1 


Vi 


1 





1 


1 





1 


1 


1 


o) 



( 

1.6779 
1.6180 
1 

1.9889 
1 

1.6180 
1 
1 

V 



1.6779 


1.6779 
1.9938 
1 

1.9843 
1 
1 

1.6532 




1.6180 
1.6779 


1.9843 
1.9889 

1 

1 

1 

1 





1 

1.9938 
1.9843 


1.7685 
1.6779 
1.6779 
1.6532 

1 





1.9889 
1 

1.9889 
1.7685 


2.4047 
1 

1.6330 
1.6482 




1 

1.9843 

1 

1.6779 
2.4047 


1.6180 
1 
1 




1.6180 
1 
1 

1.6779 
1 

1.6180 

1 
1 




1 
1 
1 

1.6532 
1.6330 

1 

1 



1.0515 




1 

1.6532 
1 
1 

1.6482 
1 
1 

1.0515 





°\ 
















0/ 



R 



' 




1.6779 


-1.4175 
0.7801 



-0.0476 
0.9951 


1.7197 
0.9989 
-0.9435 
-0.0364 
-0.3295 
-0.7671 
1.3211 
0.5331 
-0.4617 
0.8390 
-0.2880 
-0.7459 
0.3226 
0.5827 
0.1090 
0.8390 
\ 0.5331 / 



R 




2nd Moment 





a 


Special Properties 


8.89785.R 2 




1 


chiral 



312 



NATALIE ARKUS, VINOTHAN N. MANOHARAN AND MICHAEL P. BRENNER 



Packing 173 (Graph 473525): 
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*Packing 174 (Graph 479299): 
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*Packing 175 (Graph 493558): 
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*Packing 176 (Graph 506514): 
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Packing 177 (Graph 506611): 
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Packing 178 (Graph 508571): 
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Packing 179 (Graph 508895): 
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Packing 180 (Graph 508946): 
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Packing 181 (Graph 509963): 
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Packing 182 (Graph 510336): 
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Packing 183 (Graph 510412): 
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Packing 184 (Graph 512851): 
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Packing 185 (Graph 514295): 
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*Packing 186 (Graph 514375): 
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Packing 187 (Graph 516148): 
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Packing 188 (Graph 516292): 
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Packing 189 (Graph 516323): 
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Packing 190 (Graph 516437): 
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Packing 191 (Graph 520931): 
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Packing 192 (Graph 521137): 
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Packing 193 (Graph 521334): 
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*Packing 194 (Graph 550941): 
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Packing 195 (Graph 590324): 
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Packing 196 (Graph 593645): 
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Packing 197 (Graph 593702): 
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Packing 198 (Graph 593729): 
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Packing 199 (Graph 593782): 
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Packing 200 (Graph 593834): 
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Packing 201 (Graph 593842): 



A : 



V : 



/0 





1 


1 


1 





1 






















1 





1 


1 





1 























1 





1 


1 


1 














1 








1 


1 








1 














1 














1 





1 














1 


1 


1 











1 


1 














1 1 





1 











1 


1 














1 


1 


1 











1 


1 














1 


1 





1 


1 


1 





1 














Vo 1 


1 


1 1 


1 


1 


1 


1 


o) 














/ 


1 


9149 


1.7321 




1 


2.3805 


1 


1 


1.9149 


1 




1.9149 







1.7321 


1.6330 




1 


1.9149 


1 


1 


1.4142 





1.7321 


1 


7321 







1.7321 


1.4142 


1 


1.7321 


1 


1 





1 


1 


6330 


1.7321 







1.9149 


1 


1 


1.9149 


1.4142 





2.3805 




1 


1.4142 


1.9149 







1.9149 


1.7321 


1 


1.7321 





1 


1 


9149 


1 






1 


1.9149 





1.4142 


1.6330 


1 





1 




1 


1.7321 




1 


1.7321 


1.4142 





1.4142 


1 





1.9149 




1 


1 




1.9149 




1 


1.6330 


1.4142 





1 





1 


1 


4142 


1 




1.4142 


1.7321 


1 


1 


1 








V 


































0/ 



R 



' 




1.9149 


1.4434 
0.9574 


-0.0000 
0.5222 
0.8528 
0.8660 
2.1760 
0.4264 
0.8660 
0.2611 
0.4264 

-0.2887 
0.9574 

-0.0000 
0.8660 
1.6537 

-0.4264 
0.5774 
0.6963 

-0.4264 
0.5774 
1.2185 
V 0.4264 / 



R 




2nd Moment 





a 


Special Properties 


8.7000CLR 2 




1 


chiral 



DERIVING FINITE SPHERE PACKINGS 



Packing 202 (Graph 593847): 
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Packing 203 (Graph 594036): 
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*Packing 204 (Graph 594112): 
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Packing 205 (Graph 594189): 
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Packing 206 (Graph 594241): 
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*Packing 207 (Graph 594872): 
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Packing 208 (Graph 601365): 
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*Packing 209 (Graph 606541): 
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Packing 210 (Graph 607008): 
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Packing 211 (Graph 610105): 
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Packing 212 (Graph 610977): 
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Packing 213 (Graph 613779): 
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Packing 214 (Graph 618645): 
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Packing 215 (Graph 618718): 
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Packing 216 (Graph 618763): 
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Packing 217 (Graph 618802): 
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Packing 218 (Graph 618803): 
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**Packing 219 (Graph 618804): 
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Packing 220 (Graph 619044): 
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Packing 221 (Graph 619047): 
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Packing 222 (Graph 619264): 
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*Packing 223 (Graph 622056): 
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Packing 224 (Graph 622541): 
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Packing 225 (Graph 622587): 



A : 



V : 



f° 








1 


1 





1 


1 



















1 


1 





1 


1 


1 

















1 


1 





1 


1 


1 

















1 


1 


1 





1 


1 

















1 





1 





1 


1 

















1 


1 


1 





1 


1 














1 


1 


1 


1 





1 


1 




















1 


1 


1 





1 


1 








1 


Vi 


1 


1 





1 


1 


1 





1 


o) 



( 

1.4325 
1.7317 

1 

1 

1.9053 
1 
1 

1.3951 
V 



1.4325 


1.6330 
1.4730 

1 

1 

1.6667 
1 
1 




1.7317 
1.6330 


1.7318 
1.9907 

1 

1 

2.1334 
1 




1 

1.4730 
1.7318 


1.6667 
1.9237 

1 

1 

1 





1 
1 

1.9907 
1.6667 


1.6668 
1.7246 
1 

1.6328 




1.9053 
1 

1 

1.9237 
1.6668 


1.6330 
1.8846 

1 





1 

1.6667 
1 
1 

1.7246 
1.6330 


1.6328 
1 




1 
1 

2.1334 
1 

1 

1.8846 
1.6328 


1.3524 




1.3951 
1 
1 
1 

1.6328 
1 
1 

1.3524 





°\ 
















0/ 



R 



- 




-1.4325 


-1.5187 
-0.8322 


-0.1604 
-0.3080 
-0.9378 

0.3807 
-0.7163 

0.5849 
-0.9577 
-1.6343 

0.2051 
-0.9349 
-0.0958 
-0.3418 

0.5744 
-0.7163 
-0.3962 
-0.7254 
-1.0466 
-0.5699 
-0.5948 
-0.7163 
V 0.3649 / 



R 




2nd Moment 





a 


Special Properties 


8.33646^ 




1 


chiral 

new seed 



DERIVING FINITE SPHERE PACKINGS 365 



*Packing 226 (Graph 627354): 
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Packing 227 (Graph 629072): 
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Packing 228 (Graph 630908): 
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*Packing 229 (Graph 631053): 
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Packing 230 (Graph 632638): 
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Packing 231 (Graph 632698): 
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Packing 232 (Graph 632806): 
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*Packing 233 (Graph 633829): 
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Packing 234 (Graph 663203): 
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Packing 235 (Graph 665570): 
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*Packing 236 (Graph 668570): 
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*Packing 237 (Graph 679344): 
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Packing 238 (Graph 679397): 
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Packing 239 (Graph 685049): 
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Packing 240 (Graph 714961): 
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Packing 241 (Graph 730317): 
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Packing 242 (Graph 735131): 
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Packing 243 (Graph 735240): 
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*Packing 244 (Graph 736264): 
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Packing 245 (Graph 736370): 
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Packing 246 (Graph 736396): 
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Packing 247 (Graph 740549): 



A : 



V : 



f° 








1 


1 





1 


1 


1 



















1 


1 


1 


1 


1 

















1 





1 





1 


1 











1 





1 





1 





1 








1 














1 








1 


1 














1 





1 


1 


1 





1 











1 


1 





1 


1 


1 








1 


1 








1 


1 


1 





1 


1 





1 








1 


Vi 


1 


1 








1 





1 


1 


o) 



( 

1.4142 
1.4142 

1 

1 

1.7321 
1 
1 
1 

V 



1.4142 


1.6330 
1.7321 
1.9149 

1 

1 

1 

1 





1.4142 
1.6330 


2.3805 
2.1344 
1 

1.9149 
1 

1.9149 




1 

1.7321 
2.3805 



1 

2.4495 
1 

1.7321 
1 




1 

1.9149 
2.1344 

1 



2.3805 
1.6330 
1.9149 

1 





1.7321 
1 

1 

2.4495 
2.3805 


1.7321 

1 

1.7321 




1 
1 

1.9149 
1 

1.6330 
1.7321 



1 

1 





1 
1 
1 

1.7321 
1.9149 

1 

1 



1.4142 




1 
1 

1.9149 
1 
1 

1.7321 
1 

1.4142 





°\ 
















0/ 



R 








-1.4142 


-1.3333 
-0.4714 


1.0000 



0.5000 
0.2357 
0.8333 
-1.0000 
-1.4142 
-0.0000 
0.5000 
-0.7071 
-0.5000 
-0.5000 
-0.7071 
-0.5000 
0.5000 
-0.7071 
0.5000 
-0.5000 
-0.7071 
V 0.5000 / 



R 




2nd Moment 





a 


Special Properties 


9.38889i? 2 


c 2 


2 


chiral 



DERIVING FINITE SPHERE PACKINGS 387 



*Packing 248 (Graph 740809): 
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Packing 249 (Graph 743234): 
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Packing 250 (Graph 743263): 
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Packing 251 (Graph 744210): 
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Packing 252 (Graph 744235): 
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Packing 253 (Graph 744255): 
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Packing 254 (Graph 744305): 
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Packing 255 (Graph 744993): 
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Packing 256 (Graph 745020): 
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Packing 257 (Graph 749808): 
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Packing 258 (Graph 750016): 
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Packing 259 (Graph 750276): 
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Packing 260 (Graph 164942): 
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The other graphs that lead to the same 25 contact packing (packing 260), with their corresponding 
distance matrices and coordinates: 

The table containing the minimum of the 2nd moment, etc., as well as the pictures will be the 
same, and thus is not included in the list. 
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Graph 458233: 
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Graph 561797 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 594218: 
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Graph 609155 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 609244 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 609350: 
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Packing 261 (Graph 362694): 
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The other graphs that lead to the same 25 contact packing (packing 261), with their corresponding 
distance matrices and coordinates: 

The table containing the minimum of the 2nd moment, etc., as well as the pictures will be the 
same, and thus is not included in the list. 
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Graph 593947: 
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Graph 726649 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 726751: 
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Graph 726815 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 735230: 
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Packing 262 (*Graph 209285): 
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The other graphs that lead to the same 25 contact packing (packing 262), with their corresponding 
distance matrices and coordinates: 

The table containing the minimum of the 2nd moment, etc., as well as the pictures will be the 
same, and thus is not included in the list. 
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Graph 287036: 
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*Graph 457764: 
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^Graph 467242: 
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Graph 618805 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 622325 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 622375 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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Graph 622382 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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^Graph 622387: 
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+ Graph 665605 (note that as a 24 contact A this is a non-iterative graph, but as a 25 contact 
packing (V) it is iterative): 
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APPENDIX C: PSEUDO CODE 



Here we include a list of changes that were made in the code since our last paper (PRL, 
103, 118303). The entire code, both old and new versions of it, can be downloaded from the 
following website: http://people.seas.harvard.edu/~narkus/. We also include the algorithm for 
the triangular bipyramid rule, as applied to iterative packings, in which case equation 10 of the 
main text can be applied directly: 

I. Changes Made to Code 

1. In rare cases, for some of the triangular bipyramids, there are dihedral angles that are not 
defined. This occurs when >= 3 points lie on a line. This occurs in certain packings that 
contain octahedra. In these cases, the global consistency check that is normally performed 
and checks for the sum and differences of dihedral angles is not valid. In the old code, 
running this check when the dihedral angle was not defined caused some physical packings 
to be erroneously ruled out as globally inconsistent and thus unphysical. Now, the global 
consistency check (in the function 'global_consistency_check') first checks to see that all 
dihedral angles are well defined - if they are, the check proceeds as normal, if they're not 
then a special global consistency check for 3 points lying in a line (as was described in 
section 4.2 of the main text of the paper) is performed. Making this correction caused 2 
more packings to be realized at n = 9 and the bulk of the 39 more packings that were 
found at n — 10 to be realized. 

2. 'round_numi' and 'preci,' where T is a number, has been entered in various functions as 
an easier way to control the round-off numbers. Some of the round-off numbers have been 
changed so that less error is introduced as a result of these numbers. This only affects 
the number of packings found at n — 10, and accounts for around 2 more packings being 
found. 

In particular, 'prec2' in the function 'checkjround' has been changed from -6 to -10. 
This means that the precision to which we check that a certain quantity = is performed 
to 10~ 10 instead of 10~ 6 . 

3. All functions that check for 9 particle physical new seeds now first check to see that the 
distance between 2 particles, r^, is unknown before writing in the distance corresponding 
to the appropriate of that seed. (This prevents any already known distances being 
overwritten). In the old code, only 1 of the 9 particle new seed functions performed this 
check. For the functions that did not perform this check, occasionally an already known 
distance was erroneously overwritten and this in some cases caused a physical packing to 
be erroneously ruled out as unphysical. Correcting this error also only affects the number 
of packings found at n — 10; it also accounts for several more packings being found. 

2. Algorithm 

1. Loop through A. Each Aij = corresponds to an unknown r^. Keep track of which 
distances are unknown. 

2. For each unknown 

(a) Construct each possible ditetrahedron until one with all 9 distances in r is found. 
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Calculate using the triangular bipyramid rule. 

(A) If a solution exists, will fall out immediately. If < R, discard the 
solution as it is unphysical. 

(B) If a solution does not exist, this A is unphysical - move on to the next 
A. A solution will not exist if the triangle inequality is violated - this will 
show up as the absolute value of the argument of cos -1 being > 1, which 
is undefined (yielding NaN). 

Store all physically consistent r^. (These are locally consistent solutions, global 
consistency must still be checked), 
(b) If no ditetrahedron with known "F" exists, move on to the next j-y and repeat step 



2a 



3. If any unknown remain, repeat step 2a (This will be repeated until all are solved) 

4. Check for global consistency. 

The dihedral angle A\ is given by 



A\ 



(In calculating in equation 15 of the main text, we need not consider the latter two A\ 
solutions as cos(27r — x) — cos(a;)). 

For each potential set of solutions r^, construct all possible ditetrahedra in the structure. 
If 3 ditetrahedra where A\ is not satisfied by one of the above formulas, then the solution 
set is globally inconsistent, and should be discarded. 
For each iterative packing, all can be solved via this algorithm. 
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